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LOCAL WELL-POSEDNESS FOR THE BENARD CONVECTION 
WITHOUT SURFACE TENSION 

YUNRUI ZHENG 


Abstract. We consider the Benard convection in a three-dimensional domain 
bounded below by a fixed flatten boundary and above by a free moving surface. 
The domain is horizontally periodic. The fluid dynamics are governed by the 
Boussinesq approximation and the effect of surface tension is neglected on the free 
surface. Here we develop a local well-posedness theory for the equations of general 
case in the framework of the nonlinear energy method. 


1. Introduction 

1.1. Formulation of the problem. In this paper, we consider the Benard convec¬ 
tion in a shallow horizontal layer of a fluid heated from below evolving in a moving 
domain 

Q{t) = {y E S X M I -1 < ?/3 < r]{yi,y2,t)} . 

Here we assume that S = (LiT) x (L2T) for T = M/Z the usual 1-torus and Li, L 2 > 0 
the periodicity lengths. Assuming the Boussinesq approximation [1], we obtain the 
basic hydrodynamic equations governing Benard convection as 

dtu + u-Vu-\ - Vp = uAu + gaOey^, in 

Po 

dtO + u-VO = kA6, in 

u\t=o = uo{yi,y2,y3), |t=o= 6'o(yi,y2,?/3), 

Here, u = {ui,U 2 ,U 3 ) is the velocity field of the fluid satisfying divtt = 0, p the 
pressure, g > 0 the strength of gravity, > 0 the kinematic viscosity, a the thermal 
expansion coefficient, e^g = (0,0,1) the unit upward vector, 9 the temperature field of 
the fluid, K the thermal diffusively coefficient, and po the density at the temperature 
Tq. Notice that, we have made the shift of actual pressure phy p = p + gy^ — Patm 
with the constant atmosphere pressure Patm- 
The boundary condition is 

dtg-u' -Vg + us = 0, on {yg = p(t, yi, ya)}, 

(yj - z/D(tt))n = ggn + aHn + (t-V)at, on {ys = y(t, yi, y2)}, 
n-Ve + BW = -1, on {ya = y(t,yi,y2)}, 

1^3 — — 1 tl; 9 1^3 = — 1 0, 

Here, u' = (ui,tt2), I the 3x3 identity matrix, 0{u)ij = diUj + djUi the sym¬ 
metric gradient of tt, JV the upward normal vector of the free surface {ys = r/}, 
n = jVI\jV\ the unit upward normal vector of the free surface {ya = g} where 
jV = (—dig, —d 2 g, 1) is the upward normal vector of the free surface {ya = g} and 
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|.yK| = (9ir/)2 + { 821 ])^ + 1, t the unit tangential vector of the free surface, Bi > 0 

the Biot number and H the mean curvature of the free surface. For simplicity, we 
only consider the case without surface tension in this paper, i.e. a = 0 . 

We will always assume the natural condition that there exists a positive number 
5o such that 1 + 7?o ^ > 0 on S, which means that the initial free surface is strictly 

separated from the bottom. And without loss of generality, we may assume that 
po = fj, = K = a = g = Bi = l. That is, we will consider the equations 


' dtu + u ■'Vu + Vp — Au — Ocy^ = 0 

div u = 0 
dt9 + u-V9-A9 = 0 
{pi — Du)n = pn 
V0-n + 0 = -l 
u = 0, 9 = 0 

U |t=0= Uq, 9 |i=o = ^0 
dtrj + uidip + U2d2ri2 = U3 
V |t=o = Vo 


in n(f), 
in n(f), 
in n(f), 
on {ys = v{t,yi,y 2 )}, 
on {y 3 = v{t,yi,y 2 )}, (1.1) 

on {y 3 = -1}, 
in 11(0), 
on {ys = v{t,yi,y 2 )}, 
on {y 3 = v{t,yi,y 2 )}- 


The discussion of fourth equation in (ll.ip may be found in [T3]. The eighth equation 
in (|l.ip implies that the free surface is advected with the fluid. 


1.2. Previous results. Traditionally, the Benard convection problem has been stud¬ 
ied in fixed upper boundary and in free boundary surface with surface tension. 

For the problem with surface tension case, the existence and decay of global in 
time solutions of Benard convection problem with free boundary surface was proved 
by T. lohara, T. Nishida and Y. Teramoto in spaces. T. lohara proved this in 
2-D setting. T. Nishida and Y. Teramoto proved this in 3-D background. They all 
utilized the framework of [3] in the Lagrangian coordinates. 

1.3. Geometrical formulation. In the absence of surface tension effect, we will 

solve this problem in Eulerian coordinates. First, we straighten the time dependent 
domain Il{t) to a time independent domain D. The idea was introduced by J. T. Beale 
in section 5 of [3] . And in and [8], Y. Guo and I. Tice proved the local and global 

existence results for the incompressible Navier-Stokes equations with a deformable 
surface using this idea. In m, m and [8], Guo and Tice assume that the surface 
function y in some norms is small, which means 7 / is a small perturbation for the plane 
{ 2/3 = 0}- In order to study the free boundary problem of the incompressible Navier- 
Stokes equations with a general surface function 7 /, L. Wu introduced the e-Poisson 
integral method in m- In this paper, we will use the flattening transformation 
method introduced by L. Wu. We define fj^ by 

fjS _ parametrized harmonic extension ofy- 

The dehnition of can be seen in the section 1.3.1 of m for the periodic case. 
We introduce the mapping <I>'^ from D to D(t) as 


: {xi,X2,X3) {xi,X2,X3 + {I + X3)y^) = (7/1, ^2, 2 / 3 ), 


( 1 . 2 ) 
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and its Jacobian matrix 

/ 1 0 0 \ 

= 0 1 0 

\ ) 

and the transform matrix 

/ 1 0 \ 

= ((V^>^)"^)^ =01 -B^K^ 

\ 0 0 ) 

where 


A^ = (1 + X3)dir/^, B^ = (1 + X3)d2rj^, = 1 + 77 ^ + (1 + X3)5377^, = l/J^. (1.3) 

According to Theorem 2.7 in m and the assumption that l+r/o > Jo > 0, there exists 
a J > 0 such that J^(0) > J > 0 for a sufficiently small e depending on ||?7o||j|/5/2. This 
implies that $^(0) is a homomorphism. Furthermore, ‘^''(0) is a diffeomorphism 
deduced from Lemma 2.5 and 2.6 in m- For simplicity, in the following, we just 
write fj instead of rj^, while the same fashion applies to £/, $, A, B, J and K. Then, 
we define some transformed operators. The differential operators V^, div^ and 
are defined as follows. 


diVi/ u = £/ijdjUi, 

A^/ = • V^/. 

The symmetric jzZ-gradient is dehned as {0^u)ij = s^ikdkUj + s^jkOj-Ui. And 
we write the stress tensor as Ss^{p,u) = pi — where I is the 3x3 identity 

matrix. Then we note that divj^ Sj^{p, u) = for vector fields satisfying 

divert = 0. We have also written = {—dir], — 82 ^, 1) for the nonunit normal to 
{ 2/3 = 'n{yiA 2 ,t)}- Then the original equations (11.11) becomes 


' dtu - 5 * 77(1 + X3)Kd3U + u ■ V^u - A^u + V^p - = 0 in 11 

■ u = 0 in U 

dt9 — dtfi{l + X3)Kd36 + u ■ — A^9 = 0 in U 

{pi — = r/cXF on S 

< V^9-A^+ 9\A^\ = -\A^\ on S (1.4) 

u = 0, 9 = 0 on Bfo 

u(x,0)=uo, 9{x,0) = 9 q in U 

dtp + ui5i77 + U2d2P = U3 on S 

^p{x',0) = po{x') on S 
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where 63 = (0,0,1) and we can split the equation (|1.4I) into a equation governing 
Benard convection and a transport equation, i.e. 


dtu - (9t7y(l + X 3 )Kd 3 U + u ■ - A^u + V^p - 9V^y^ = 0 

in 

n 

• ti = 0 

in 

n 

dtO - dtffil + x 3 )Kd 3 e + u ■ - A^e = o 

in 

n 

{pi — = rf^yV 

on 

S (1.5) 

+ e\jT\ = -|./C| 

on 

s 

II 

II 

0 

on 


ti(a;,0) = uq, 9(x, 0) = 9o 

in 

n 


and 

{ dtrj + uidiTj + U 2 d 2 ri = U 3 on S 
rj{x', 0 ) = r]Q{x') on S 

Clearly, all the quantities in these two above systems are related to 77 . 


1.4. Main theorem. The main result of this paper is the local well-posedness of the 
Benard convection. Before stating our result, we need to mention the issue of com¬ 
patibility conditions for the initial data {uo,0o,r]o). We will study for the regularity 
up to N temporal derivatives for N > 2 an integer. This requires us to use uq, 6 q 
and ryo to construct the initial data 5jtt(0), d( 6 { 0 ) and d^r]{ 0 ) for j = 1 ,..., and 
5^p(0) for j = 0,..., N — 1. These data must then satisfy various conditions, which 
we describe in detail in Section 5.1, so we will not state them here. 

Now for stating our result, we need to explain the notation for spaces and norms. 
When we write \\dlu\\jjk, and \\dlp\\jjk, we always mean that the space is 

and when we write ||c^r 7 ||//s, we always mean that the space is where 

and are usual Sobolev spaces for fc, s > 0. 

Theorem 1.1. Let N >2 he an integer. Assume that ryo + 1 ^ > 0, and that the 

initial data {uo, 6 o,r]o) satisfies 


■— ||'^^o||^ 2 iv + || 0 o||^ 2 jv + ||%||^ 2 iv+i /2 < 00 , 


as well as the N-th compatibility conditions (|5.28p . Then there exists a 0 < Tq < 1 
such that for any 0 < T < Tq, there exists a solution {u,p, 6 ,r]) to (|1.4n on the interval 
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[0, T] that achieves the initial data. The solution obeys the estimate 


N 


^2iV-2j + \\d{u\\‘^2fj2N-2j + l 


,=0 \0<*<^ 
N-1 / 


+ II^PII^2JV-2j-1 + \\d^p\\‘j^ 2 H 2 N- 2 j 


,=0 \0<*<^ / 

+ Z] f ®^P \\9t^\\H^N-2j + \\d{d\\l2j^2N -2j + l ]+\\dr^0k^^y 

f^Q\0<t<T ) 

( ^ 

+ sup ||??||^2iV+l/2(s) +yi sup \\dl7^\\]j2N-2j+i/2 


0<t<T 
N+1 


0<t<T 


+ \\dtP\\‘l^2fj2N-2j+S/2 

3=2 

<C(Oo,<5)P(^o), 


(1.7) 


where C{Qo,6) > 0 depends on the initial domain Hq o.nd 5, P{-) is a polynomial 
satisfying P{0) = 0, and the temporal norm Lf is computed on [0, T]. The solution is 
unique among functions that achieve the initial data and for which the left-hand side 
of (II.7p is finite. Moreover, i] is such that the mapping ^{-p) defined by (II.2p is a 
(J2N-1 jQTf each t S [0, r]. 

Remark 1.2. The space is defined in section 2 o/ [6]. 

Remark 1.3. Since the mapping is a diffeomorphism, we may change 

coordinates to produce solutions to m- 

1.5. Notation and terminology. Now, we mention some definitions, notation and 
conventions that we will use throughout this paper. 

1. Constants. The constant C > 0 will denote a universal constant that only depend 
on the parameters of the problem, N and 11, but does not depend on the data, etc. 
They are allowed to change from line to line. We will write C = C{z) to indicate 
that the constant C depends on z. And we will write a <b to mean that a < Cb 
for a universal constant C > 0. 


2. Polynomials. We will write P{-) to denote polynomials in one variable and they 
may change from one inequality or equality to another. 

3. Norms. We will write for P[^{Q) for k > 0, and R*(S) with s G M for usual 
Sobolev spaces. Typically, we will write = L^, With the exception to this is 
we will use L^([0,T];R^) (or L^([0, T]; R^(S))) to denote the space of temporal 
square-integrable functions with values in (or 77®(S)). 

Sometimes we will write || • ||fc instead of || • or || • We assume that 

functions have natural spaces. For example, the functions u, p, 6 and fj live on n, 
while T] lives on S. So we may write || • \\jjk for the norms of u, p, 0 and f] in n, 
and II • II//S for norms of rj on S. 
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1.6. Plan of the paper. In section 2, we develop the machinery of time-dependent 
function spaces based on [6]. In section 3, we make some elliptic estimates for the 
linear steady equations of (|1.8I) . In section 4, we will study the local existence theory 
of the following linear problem for {u,p,6), where we think of p (and hence £/, 
etc.) is given: 

' dtu - - 6»V^?/3 = 

•n = 0 

dte - A^e = 

(pi = F^ 

+ e\jK\ = F^ 
u = o, e = o 


in Q, 
in fl, 
in fl, 
on S, 
on S, 
on Sfe, 


( 1 . 8 ) 


subject to the initial condition tt(0) = uq and 0(0) = 9q, with the time-dependent 
Galerkin method. In section 5, we construct the initial data and do some estimates 
for the forcing terms. In section 6, we construct solutions to (|1.4p using iteration and 
contraction, and complete the proof of Theorem 11.11 


2. Functional setting 

2.1. Function spaces. Throughout this paper, we utilize the functional spaces de¬ 
fined by Guo and Tice in section 2 of [6]. The only modification is the definition of 
space For the vector-valued space its definition is the same as [6]. 

The following is the definition for the scalar-valued space 

^\t) :={0|||0||^i <oo,0|s, = O} 

1 /2 

with the norm \\0\\^i := (0,0)where the inner product (•, •)is defined as 

(0,(/))^i := [ (V^(i)0 • V^(i)0) J(t). 

Jn 

The following lemma implies that this space is equivalent to the usual Sobolev 
space H^. 

Lemma 2.1. Suppose that 0 < eo < 1 \\t] — < ^o- Then it holds that 

l|9|ll. + II>IoIIhw«(e)) ll«lll.. (2.1) 

3 S X S (i + lkollH.«,E))‘‘ l|e|ll.(n,. (2.2) 

Proof. From the Poincare inequality, we know that ||0||//i is equivalent to ||V0||j:^o. 
So in the following, we will use ||0 ||hi instead of ||V0||j:^o. 

From the assumption and the Sobolev inequalities, we may derive that 

d ^ ^ 1 + II0||l°° + ^ 1 + II0|Ih 5/2 ^ 1 + ll0ol|jy5/2. 



< max{l, \\AK\\l^,\\BK\\l^,\\K\\l^} 

< 1 + (1 + \\Vp\\l^)\\K\\l^ < (1 + \\po\\H^/2f . 


and 
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Thus (12.111 is clearly derived from the estimate of || and we have that 


■/|V^0p<(l + ||r?o||/^5/2) / 


\BK\ 


K\ 




^ (1 + ll%lli75/2) max{l, \\AK\\l, 

Now we have proved the second inequality of (12.2p . 

To prove the first inequality of (12.211 . we rewrite the \\0\\^i as 

[ J|V^0|2= [ J|V^„0|2+ [ j(v^0 + v^,0)-(V^0-V^3 0), 
J Q, J Q, J ^ 


Here £/q is in terms of r/o- By the estimates of 


/ 

Jn 




>_ 

~ 1 + 


1 


775/2 Jq 


we derive that 


1 + ll%ll775/2 7^0 


|V(0ocl.(O))|^ 


>_ 

~ (1 + 


1 


|J70 


1/75/2)^ 

where in the last inequality, we have used the following Lemma l3.ll since <h(0) is a 
diffeomorphism. Here Jq is in terms of r]Q. Then, using the estimates of ||j 2 ^'||j;,oo and 
1 1 , 00 , we have that 


/ 


J{V^6 + • (V^0 - V^,e) 


< 




(1 + ||%||j 75/2) Il^l|j7i- 


m||L° 


|J71 


Then taking eq sufficiently small, we may derive that 

^ f f J(y^0 + V^,0)-iV^0-V^,0) 


'n 


> _ 

~ (1 + 

This is the first inequality of 


1 


775 / 2 )" 


|771- 


□ 


We define an operator ICt by ICt0 = K{t)0, where K{t) := K is defined as (|1.3I1 . 
Clearly, ICt is invertible and IC^^Q = K{t)~^Q = J{t)Q, and J{t) := J = 1/K. 

Proposition 2.2. For each t G [0,T], ICt is a bounded linear isomorphism: from 
H^{Q) to H^{Cl) for A; = 0,1, 2; from L'^{Cl) to J^^{t); and from oH^{d) to J^^{t). 
In each case, the norms of the operators ICt, ICf^ are bounded by a polynomial 
P(||r/(t)||^7). The mapping 1C defined by IC0{t) := ICt0{t) is a bounded linear isomor¬ 
phism: fromL‘^{[Q,T]-H’^{d)) to L‘^{f),T];H^{d)) fork = 0,1,2;/ram ^^([0, T]; 
to and from oH^{Cl) to . In each case, the operators 1C and IC~^ are bounded 
by the polynomial P(supo<f<T Mt)\\^r). 

Proof It is easy to see that for each t G [0,T], 

11^7^11770 <||/C/|bo||0||77O<B(||,7(t)|| 


7 

777 


770 ) 


(2.3) 
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ll^i^ll//2 < ||/Ct||cl ||0||_H-2 + ||/Ct||j^2||0||c-O < P{\\'>l{t)\\^7)\\0\\H^. (2.5) 

These inequalities imply that ICt is a bounded operator from to , for fc = 0,1, 2. 
Since K,t is invertible, we may have the estimate ||/Cr^0||j|^i: < P{\\'ri{t)\\ ? 

Thus, Kt is an isomorphism of to H^, for k = 0,1,2. With this fact in hand. 
Lemma l2.ll implies that JCt is an isomorphism of L^(n) to and of oP[^{Q,) to 

The mapping properties of the operator JC on space-time functions may be estab¬ 
lished in a similar manner. □ 


2.2. Pressure as a Lagrange multiplier. The introduction of pressure function 
has been studied by Guo and Tice in section 2 of [6], of which the modification was 
given by L. Wu in section 2.2 of m- So we omit the details here. 


3. Elliptic estimates 

3.1. Preliminary. Before studying the linear problem (|1.8I) . we need some elliptic 
estimates. In order to study the elliptic problem, we may transform the equations on 
the domain H into constant coefficient equations on the domain fl' = where $ is 

defined by (|1.2I) . The following lemma shows that the mapping <1> is an isomorphism 
between P[^{Q') and Here, the Sobolev spaces are either vector-valued or 

scalar-valued. 

Lemma 3.1. Let T : H —)• 11' 6e a (7^ diffeomorphism satisfying T S , VT —/ S 
and the Jacobi J = det(VT) > 5 > 0 almost everywhere in LI for an integer 
k > 3. If V £ then u o T G for m = 0, 1,... , A; -|- 1, and 

||u o < C (^IIVT - I\\H>‘{n)^ 

where GdlVT — is a constant depending on ||VT — I\\f{k(^Qy Similarly, for 

u G we have u o G for m = 0,1,... ,k + 1, and 

\\u o < C (^||VT - \\u\\H"^{n)- 

Let S' = 'I'(S) be the top boundary ofLl'. If v £ iA™'“^(S') for m = 1,... ,k — l, then 
uoT G iL™"3(S), and 

If u £ (S) for m = 1,... ,k — 1, then u o G i /™“2 (S') and 

Ik o ~ ^ (ll^^ “ lkllj|/"*-7(S)’ 


Proof. The proof of this lemma is the same as Lemma 3.1 in [6], which has been 
proved by Y. Guo and 1. Tice, so we omit the details here. □ 
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3.2. The j 2 /-stationary convection problem. In this section, we consider the 
stationary equations 


div^ F^(p, u) - 6'V^y3 = F^ 

in 

PL 

div^ u = F"^ 

in 

P 

-A^9 = F^ 

in 

P 

S^{p,u)PP = F^ 

on 

S 

V^9-jP + 9\J^\=F^ 

on 

S 


u = 0, 9 = 0 on Sfe 


(3.1) 


Before discussing the regularity for strong solution to (13.11) . we need to define the 
weak solution of equation (13.11) . Suppose G G , F^ G 

G F“ 2 (S) and F^ G H~ 2 {T,), {u,p, 9) is called a weak solution of equation (|3.1I) 
if it satisfies • u = 


and 

^ (D^u,B^V')^o + (p, - (6'V^y3,V’)^o = (F\V')(jri)* - ’ 

(3.3) 

for any (^, '0 £ 


Lemma 3.2. Suppose F^ G F^ G F^ G (^i)*, F^ G F"i(S) and 

F^ G F~ 2 (S). Then there exists a unique weak solution {u,p,9) G x x 
to (1311) . 


Proof. For the Hilbert space with the inner product (0, 0) = (V^0, V^(f )+ 
(0 IcxFl ,0)j:^o(s)) we can define a linear functional G by 

for all 4> £ Then by using the Riesz representation theorem, there exists a unique 
9 G such that 


(V^6(, V^0)j^o + {9 \ jV\ , 0)j^o(s) — (-^^j ’ 


for all 0 G 

By Lemma 2.6 in [^, there exists a u G such that div^ u = F^. Then, we may 
restrict our test function to 0 G . A straight application of Riesz representation the¬ 
orem to the Hilbert space with inner product defined as (u, 0) = (B^u, B^0) 
provides a unique w G such that 


^ (B^u;,B^0)^o = (B^u,B^0)^o + (6'V^y3,0)^o + (F\0)^^i^* 

for all 0 G . Then we can find u satisfying 




(3.4) 


\ (B^u,B^0)^o - (0V^y3,0)^o = (i^\0>(^i)* - ’ 


by u = to + u G with div^ u = F2. 


(3.5) 
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It is easily to be seen that u is unique. Suppose that there exists another u still 
satisfies (j3.5l) . Then we have div^(ii — 11) = 0, and (]D)^(ii — u),= 0 for 
any tp ^ . By taking pj = u — u, and using the Korn’s inequality, we know that 

||rt — ■uIIji^o = 0 which implies u = u. 

In order to introduce the pressure p, we can dehne A G as the difference of 

the left and right hand sides of (|3.4p . Then \{'tp) = 0 for all ip G ^. According to the 
Proposition 2.12 in [13], there exists a unique p G satisfying (p, div^ 'ip )= X{'ip) 
for all Ip G □ 


In the next result, we establish the strong solutions of (13T]) and present some 
elliptic estimates. 


Lemma 3.3. Suppose that rj G (S) for k > 3 such that the mapping <I> defined 

in (11.211 is a diffeomorphism of kl to Q' = <I>(II). If G G G , 

F^ £ H 2 and F^ G then the problem (|3.1[) admits a unique strong solution 

{u,p,9) G X H^{Q) X i.e. {u,p,6) satisfy (|3.1[) a.e. in 12 and on S, 

Sfe. Moreover, for r = 2,... ,k — I, we have the estimate 


\u\\h^ + IIpIIh'’'-! + II^IIh’' ^ C{r])(\\F^\\^r-2 + ||F^||j^r-l + ||F^||j^r-2 


+ IIF II 3 +IIF'"II 3 


(3.6) 


whenever the right-hand side is finite, where C{if} is a constant depending on 




Proof. First, we consider the problem 


-A^9 = F^ in n, 

< V^9-M^ + 9\M^\ =F^ on S, 
0 = 0 on Sft. 


Since the coefficients of this equation are not constants, We transform this problem 
to one on 12' = <h(f2) by introducing the unknowns 0 according to 0 = 0o$. Then 0 
should be solutions to the usual problem on 12 ' = {—1 < < ??(yi,?/ 2 )} with upper 

boundary S' = {y^ = 77 }: 


' -A0 = F^o $-1 = in 12', 

< V0-^ + 0|=yF| =F^o$-i = on S', 

, 0 = 0 on Sj. 


(3.7) 


Note that, according to Lemma [3Tl G and G^ G 27^/^(S'). Then we may 

argue as the Lemma 2.8 in | 2 ] and use the Theorem 10.5 in [T], to obtain that there 
exists a unique 0 G 77^(12'), solving problem (|3.7I) with 

Il©k2(f20 < CmiG^HOin') + 

for C{r]) a constant depending on ||r 7 ||^j,+ i. 
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(3.8) 


For the .e/-Stokes equations, we introduce the unknowns v,q hy u = u o and 
q = po For the usual Stokes problem 

' S{q, v) - 063 = F^o ^>-1 = in O' 

V • u = o in O' 

S{q, v)FK = F^o = G* on S' 
u = 0 on Sfe, 

we use the same argument as in the proof of Lemma 3.6 in [6] with G^ + 063 instead 
of G^. Then we have that there exist unique v G F^(O'), q G F^(O'), solving problem 
()3.8p with 

Il^ll/i2(n/) + ||g||/ii(n') ^ C{ri) (^||G'^||/^o(q/) + IIG^IlHi(n') + + ll®llH0(n') 

for G{p) a constant depending on 




. so we have that 


+ Ikllfl-qo') + I|0 ||h2(o/) < C'(?/)(^||G'^ 11/10(11') + \\G'^\\m{n') 


(3.9) 


^fc+i • Then we may argue it as in Lemma 3.6 


for G{p) a constant depending on 
of [6] to derive that, for r = 2,. .., A; — 1 

lbll//'-(fi') + Ikll//’--i(fi') + l|0||//'■(!^^') 

^ C'(^)(||G'^II//’—2(0') + l|G'^ll//'-i(o') + I|G'^II//'-2 (o') (3.10) 

for G(r/) a constant depending on ||/?||^fc+i. 

Now, we transform back to fl with u = uo<h,p = go<h and 0 = 0 o <h. It is readily 
verified that (u,p, T) are strong solutions of (13.ip . According to Lemma [3Tl 

+ IIpII//’—1 + ^ G{p)(\\F^\\ur-2 + ||F^||j^t—1 + ||F^||j^t—2 




+ F 3 + F°ll 3 

" "//’'-i(s) " "//’-i(s)F 

whenever the right-hand side is finite, where G(ri) is a constant depending on ||r?|| 

This is what we want. □ 


In the next lemma, we verify that the constant in (|3.6I) can actually only depend 
on the initial free surface. 

Lemma 3.4. Let k > 3 be an integer and suppose that p G H ^'^2 (S) and po G 
iL^+ 2 (S). Then there exists a positive number eq < 1 such that «/||i? —?7o|| k-^ ^ 
the solution to dST]) satisfies 


miH^ + mlHr-i + 


\w 


< 


C{P0)(\\F^\\Hr-2 -I- ||F^||j^r-l -I- ||F^||j^r-2 


F^ 3 F° 3 

" "//’"-i(s) " "//’-i(E);’ 


(3.11) 


for r = 2,... ,k — 1, whenever the right hand side is finite, where G{po) is a constant 
depending on ||?7o||^fc+i . 
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Proof. Here, we use the same idea as in Lemma 2.17 of m- We rewrite the equation 
()3.1I) with its coefficients determined by r/o, i.e. it can be thought as a perturbation 
of equations of (13.11) in terms of initial data, 

div^o u) - 2/3,0 = 

V^,-u = F^ + 

-A^^e = F^ + 

S^,{p,u)pKo = F^ + F^^^ 

V^,9-pK> + 0\^o\=F^ + F^’^ 

u = 0, 9 = 0 

where 

F^'° = • Ss/{p,u) + • Ss^Q-^{p,u) + 6 (V^o _^?/3 + 6 (V^o(2/3,0 - 2/3), 

F^’° = div^o_^ n, 

F^’'^ = 

F‘^’° = F^o(p,u)(^ - JP) + 

• (=yKo - ^) + • ^ + 0 (l^ol - 1^1) • 

Here, 32 ^ 0 , =ylo and 1 / 3,0 are quantities of 32 /, jV and 2/3 in terms of t/q. By the assump¬ 
tion, we know that r/ - r/o G and \\r] - r]of ^_3 < \\p - po\\ .s < 1 

for any positive integer i. By the straightforward computation, we may derive that 

Wf^'^Wh’^-^ < C" (^1 + Iholl^fc+i) Wv - voW^k-^ dkllH’- + IIpIIj/’—1 + II^IIh-’-- 2 ) > 

\\F‘^’^\\Hr-i < c(l + ||r/o||^fc+i ) \\t] - %||^fc-3 \\u\\h-, 

\\F^’^\\Hr -2 < C(l + llr/oll^fc+i ) \\r] - mW^k-^\\0 \\h-, 

||F^’°||^^_3^^^ < C(^l + ||r/oll^fc+i ) h - 9o\\^k-3 {\\u\\h- + WpWh-^) , 

+ ll^oll^.+ i )'h - \\9\\Hr, 

for r = 2,..., A; — 1. 

Based on the Lemma 13.31 we have the estimate 


in H 
in H 
in H 
on S 
on S 
on Sfe 


(3.12) 


kll/i''- + IIpII_h^-i + II^IIh"^ 

< Cirjo) (||Fi + F1’0||^.-2 + \\F^ + F'^^^Wjjr-i + jjF^ + F 


'3,01 


\H^-2 


+ ||F^ + F' 


4,01 


3 + Wpi’ _L p 


5,01 
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where C{t]o) is a constant depending on ||? 7 o|| 4 .+ 1 . Combining the above estimates, 
we have 

< Cirio)(^\\F^\\Hr-2 + \\F^\\h-^ + \\F^\\m-^ + 

+ C'(r/o)(l + ||r/o||^fe+i) \\v-Vo\\^k-^ {\W\\Hr + \\p\\h-^ + Mh-) , 

(3.13) 

for r = 2,... ,k — 1. Then, if \\r] — ryoll^fe-l is to be chosen small enough such that 

the second term of the above inequality on the right-hand side less than |(||u||//r -|- 
1 + II^IIh^)) then it can be absorbed into the left hand side, and we have that 


+ IIpIIh’--1 + \W\\h^ 

< Cim){\\F^\\Hr-2 + \\F^Hr-l + \\F^\\Hr-2 + 11 I + 11 | 




for r = 2 ,..., /c — 1 . 


□ 


Notice that the estimate in ()3.1ip can only go up to /c — 1 order, which does not 
satisfy our requirement. In the next result, we can achieve two more order with a 
bootstrap argument, where we use the idea of m- 


Proposition 3.5. Let k > 3 be an integer. Suppose that p € as well as 


po S (S) satisfying \\p - po\ 


H 


< Co- Then the solution to (13.ip satisfies 

^ (S) 


WWh^ + IIpII 




-1 + 




< 

rs^ 


C{po){\\F^\ 


H 


-2 + IIF^ 


\H 


■-1 -|- IIP II rfr-2 + IIP 


I 3 + IIP II 3 


(3.14) 

for r = 2,..., k + 1, whenever the right hand side is finite, where C{po) is a constant 
depending on ||?7o||^fc+i 


Proof. Here, we only consider the case for r = k and r = k + 1, since the conclusion 
has been proved when r < k — 1. For m G N, we define p'^ by throwing away high 
frequencies: 


p"^{n) 


p{n), for |n| < m — 1 , 
0 , for |n| > m. 


Then for each m, r/™ G H^{YL) for arbitrary j > 0 and p"^ —)• 77 in as 

m —>■ 00. 

We consider the problem (13.ip with 32 / and ,JV replaced by and , and 773 
replaced by y™. Since p'^ G if 2 , we may apply Lemma 13.31 to deduce that there 
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exists a unique 0™) which solves 


ly™ = in 0. 


div^m u™ = in 
-A^mO^ = F^ in n 


= on E 
^ QTU \^m\ ^ F^ On S 
u^ = 0, 0^ = 0 on Eb 


(3.15) 


and satisfies 
Wu^lnr + Up" 


Ij^r-l 


+ 11^ 


mil cT' 

Ur ^ 

iJ /-v^ 


C(||p" 


i]jk+^ 


5)(||f^II^.-2 + iif^ 


I//’’ 


-1 + 11-^ 


+ llF'‘ll 3 +llF^II 3 

for r = 2,..., A: + 1. In the following, we will prove that the constant C'(||ry™'||^j,^5) 


can be improved only in terms of 
For convenience, we define 


'h^+\ ' 


!^ = C{r^o)P{0^)(\\ 


7?1||2 _l II 772112 I II 773 1|2 i || p4||2 

.C \\jj, — 2 + ||.r Wfj, — 1 + ||n \\fjr-2 F W-T 


\Hr-i 


\Hr-2 




+ \\FY 3 


where C'(yo) is a constant depending on 


H 


k+i 


and P{t]) is a polynomial of 


Then after the same computation as in the proof of Proposition 2.18 
in m except for the only modification of F replaced by F^ + we have 




for r = 2,... ,k + 1. That’s because in the above estimate, we only need to consider 
the terms || 0 ”^||/ 7 '’’, for r = 2, • • • ,k — 1, but || 0 ”^||/ 7 '’’ < is assured by the Lemma 

ETH 

Then we consider the temperature 0"*. In the following of bootstrap argument, we 
may abuse the notation 9 instead of 0 ™ and also for 77 , £/, but they should be 
thought as y™, j?/™, c/F™. We write explicitly the equation of 9 as 

5ii6» + 8220 + (1 + 742 + B'^)K‘^d339 - 2AKdi39 - 2BKd239 

(3.16) 

+ {AKd3{AK) + BKd3{BK) - di{AK) - d2{BK) + Kd3K)d39 = -F^. 


step 1 r = k case. By Lemma 13.111 

ll«lll,.-. < C(,„) (||f3||^._, + I|F=II^._5,,,,) < S", 

where the constant C{t]q) only depends on Hyoll /t+i • * = 1)2, since di9 

satisfies the equation 

—Aj^di9 = F^ in 17, 

< V^di9 • ^ + 9*0 1^1 = F^ on E, 

9*0 = 0 on Eb, 


where 


F^ = diF^ + diva*^ V^0 + div*^ Va.^0, 

i75 = Q.p^ _ ^ _ v^0 • 9*^ - 09* 1^1. 
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Applying the Lemma A.1-A.2 in [^, we have 


\F- 


< 


3||2 


\Hk-3 


-3 + \\F 


5||2 




IF 


'3||2 


\Hk-2 


+ l|F 


5||2 




2 

Hk-i 




Employing the k — 1 order elliptic estimate, we have 

\m\m-, < C(,,„)(||C||?,„ + lir'IlL-lig) s 

Then taking derivative on both sides of (13.161) and focusing on the term 
(1 + + B‘^)K‘^d^6, the estimates of all the other terms in Lf^-norm implies 

that 

ll^t^llio < 

Thus, we have proved that 

ll^ll^. < 

step 2 r = k + 1 case. 

For i,j = 1,2, since dijO satisfies the equation 

-A^dijO = in 11, 

< V^dijO • ^ + dijO I = F^ on S, 
dijO = 0 on Sfe, 

where 


p3 

_p5 


= dijF^ + diva,^.^ V^e + div^ Vg, .^6» + div^,^ + divg^,^ Va,^6» 

+ divg.^ V^djO + diva^.^ V ^diO + div^ ^ dis^djO + div^ 

= d^jF^ - v^e • - (Va,^0 + v^d^6) • - {Vg.^e + 

- + Vg^^djO - Vg.^d^G)^ - mj |=A^| - di0dj 1^1 - djm i^|. 

Applying the Lemma A.1-A.2 in [6] to the forcing terms, we have 

|^5||2 


I ||2 I 




< 


I 7?^II2 I 

|.C \\j^k-l + 


1 




Then the Lemma l3 .11 1 implies that 


lia««llL-. < cOo) (iT'llL-. + iT'^llL-i,,,,) ;£ 

Since we have proved the case r = k, we take derivative d^~‘^di on both 
sides of (I3.16P for i = 1,2 and focus on the term of (1 + + B‘^)K^d^di9. 

Utilizing the estimates of all the other terms in Lf^-norm, we have 

wdlmiHo < 

Then, taking derivative on both sides of (j3.16p and focusing on the term 
of (1 + A^ + B‘^)K‘^d 2 ^^ 0 , by all the estimates above, we have 
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Therefore, we have proved 

\mW < 


Now, we go back to the original notation. According to the convergence of r/™, we 
have 


\u 


IH’ 




C{m)P{rin{\\F^fHr-. + \\F^ 


12 ill Z7'3 ||2 ill 77^4 II2 i m ^5 ||2 


Civo)P{v){\\F^\\lr-2 + WF^Wlr-^ + \\F^\\lr-2 + ||T"||L _3 


C(r?o)(||i"' 


\H 


l^r—2 ~r ||-t llj^r—1 “r ||-t- lljf'^ 

_L II 7 ?^ 11 ^ _L II 7 ? 3||2 I II r ^4 

^_2 + ll-T llj^r-1 + W-t^ \\ffr-2 + 11 


+ l|i^®llL^3 


IH'’- 




++"5(S) 

+ ||^5„2 


(3.17) 




for r = 2,... ,k + 1, where in the last inequality we have used the assumption that 
\\r] — ??o||„fe+i < £o and the term P(t/o) is absorbed by C{t]q). Here C{r]o) depends 
only on ||r7o||^fe+i. 

The inequality of boundedness (I3.17jl implies that the sequence {(tt—,p—, 0—)} 
is uniformly bounded in X x 77*', so we can extract a weakly convergent 

subsequence, which is still denoted by {(tt—,p—, 0—)}. That is, u"* ^ in 77^(H), 
pm pO 77^“i(17) and 0— ^ 0*^ in 77^(0). Since r/— —)■ r/ in we also 

have that j?/— —)> si., J— —)■ J in 77^(11), and —)• jV in 77^“2(S). 

After multiplying the equation div ym _ p 2 j^y uuJ^ for w G C'“(n) and 
integrating by parts, we see that 


[ F^wJ^ = [ div^^(u-)u;J- = - [ u- • 

Jq Jq Jq 


Jo, 


/ divS!^{u^)wj, 
JQ 


from which we deduce that div^ vP = F"^. Then multiplying the third equation in 
()3.ip by tcJ— for w G o77^(H) and integrating by parts, we have that 


[ v^m0- • + [ e^w |+b-| = [ 

Jq, Jt, Jq 


= / F-^wJ^F 




which, by passing to the limit m ^ oo, reveals that 


[ V^0°-V^u;J+ [ 0°u;l^|= [ . 
J Q t/ S J Q 


F'^wJ + 


F^w. 


Finally we multiply the first equation in (|3.1|) by tcJ— for w G o77^(17) and integrate 
by parts to see that 


- / F^-w. 


[ - p-J- - • wJ^ = [ PP wr 

Jq 2 Jn 

Passing to the limit m ^ oo , we deduce that 

— / : D^u;J+ p° div^(rt;) J — 0'^V^y3 • rcJ = f F^-wJ— / F^-w. 

Jq 2 Jn Je 


After integrating by parts again, we deduce that (ri°,p°,0°) satisfies (13.ip . Since 
{u,p,9) is the unique solution to (j3.ip . we have that u = u^, p = and 0 = 9^. 
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Then, according to the weak lower semicontinuity and the uniform boundedness of 
()3.17p . we have that 


+ IIpIIh'—1 + 
<Civo){\\F^ 




H '—2 + II cjr-l + ||T^||i/ 7—2 + ||T''^|| ^3 + ll-F^II 3 

^ . . . . . 


for r = 2,..., /c + 1, where C{r]Q) is a constant depending on 


J7''+2(S)' 


□ 


3.3. The j 2 ^'-Poisson problem. Now we consider the elliptic problem 


' A^p = f 

inll. 


P = f 

on S, 

(3.18) 

, V^p- z/ = f 

onSfe, 



where ly is the outward-pointing normal on The details of elliptic estimates of 
(|3.18l) has been interoreted in [U] and m, so we omit them here. 

4. Linear estimates 

Now we study the problem (11.81) . following the path of [6]. First, we will employ 
two notions of solution: weak and strong. 

4.1. The weak solution. Suppose that a smooth solution to (II.8p exists, then by 
integrating over 11 by parts, and in time from 0 to T, we see that 

{dtu,'ip) ^ (w,V’)L2jri - (P>div^iA)^2^o - (^V^ys, V’)l2^o 

+ (V^0, V^(/>)^2^o + (6*1^1 A)l'^h^[Ts) 

for (j), if; G 

If we were to restrict the test function ijj to G 3F, the term (p, div^ 
would vanish. Then we have a pressureless weak formulation. 


(4.1) 


(4.2) 


{dtu, 'ip)j^2jg>o + - {u, V’)i,2j^i — {OV^ys, V’)i,2^o 
~ “ {F )'*/’)i2jyO(S) ! 

idt0,(t))j^2j^o + (V^6(, V^(/!))^2j^o + (^l-^l A)l'^h°{-e) 

= {F^,'4^)i2^o + {F^^'^)l'2hO{T.) ’ 

This leads us to define a weak solution without pressure. 

Definition 4.1. Suppose that uq G ^{0), Oq G if°(n), G {^t)* and F^ + 

F^ G . If there exists a pair (u, 9) achieving the initial data uq, Oq and satisfies 

u G Jifr^, 6 G and dtu G , dt9 G , such that 


+ - (u,V’)l2^i - (0V^y3,V’)L2^o = [F^ - , 


(4.3) 


{dt9,4>){^f)* + (^! + {G \^\ ) 4’)l‘^h°{t.) — {F^ + F^ ,fi) , 

holds for any G and (j) G , we call the pair {u, 9) a pressureless weak solution. 
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Since our aim is to construct solutions with high regularity to (|1.8I) , we will directly 
construct strong solutions to (14.3p . And it is easy to see that weak solutions will arise 
as a byproduct of the construction of strong solutions to (II.81) . Hence, we will not 
study the existence of weak solutions. 

Now we derive some properties and uniqueness of weak solutions. 

Lemma 4.2. Suppose that u, 9 are weak solutions of (14.3p . Then, for almost every 
t G [0,T], 


^ 0 


+ S f f \u{s)fdsJ{s)ds+ [ [ 9{s)V^y3-u{s)ds, 
^ Jo Jn Jo Jn 


l|l«(«)ll3r.(,) + /ll«Wll3r.,.)d» 


+ 



+ {F^ + + - 



0 Jn 


|6'(s)|^clsJ(s)ds. 


Also, 


(4.4) 


sup \\e{t)\\%.o^t) + ll^ll^i ^ {Co{t])T) 

0<t<T T 


2 ill p3 pS||2 

M’^io) IH “T ll(jri)* I ’ 


(4.5) 


sup ||u(t)||^o(t) + \\u\\%,i < exp {CCo{v)T) (||m(0)||^o(o) + 11^(0)||^o(o) 

0<1<T 

where Co{r]) := max{supo<t<T \\dtJK\\L^,supQ<t<T l|V^2/3||i,°°}- 

Proof. The identity (|4.4p follows directly from Lemma 2.4 in [B] and (14.3p by using 
the test function tp = ux[o^t] G ^T, and (p = 9x[o,t\ £ where X[o,t] is a temporal 
indicator function to 1 on the interval [0,t]. 

From p4.4l) . we can directly derive the inequalities 

2 ll^(^)ll^o(t) + ll^ll^d - 2 “*■ 11-^^ “*■ + 2^o(i?)||6'(i)||^^o, 

(4.7) 

+ 2^o(^)ll'“(^)ll^o + ^^0(^)116*11^^^111^11^^1, 


where, for (14.8p . we have used the Poincare inequality in Lemma A. 14 on [6], and 
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and similarly for - F^||(jrt)*, 11^^ + Inequalities (I4.7|) . (|4.8I) 

and Cauchy inequality imply that 

2p{'t)\\%’0{t) + - 2 + 11 -^^ “ + 2^o{'n)P{t)\\%’0, 

2 ll'“(^)ll^o(i) + gll^^ll^^i < 2ll'“(^)ll^o(o) + 11-^^ “ + 2^o{'n)\W{'t)\\%’0^‘^'^'^ 

+cCo{v)\\ef^u 

Then (|4.5I) and (14.61) follow from the integral inequality (14.91) and Gronwall’s lemma. 

□ 

Proposition 4.3. Weak solutions to (14.3p are unique. 

Proof. Suppose that (u^, 9^) and (u^, 9“^) are both weak solutions to (14.3p . then {w, i?), 
defined hy w = — u'^ and '9 = 9^ — 9“^, is a weak solution with = 0, 

F^ + F^ = 0, ui(0) = M^(0) — tt^(O) = 0 and = 0^(0) — 0^(0). Then the bounds 
dM]) and ()4.6p imply that w = 0 and i? = 0. Hence, weak solutions to (14.3p are 
unique. □ 

4.2. The strong solntion. Before we define the strong solution, we need to define 
an operator Dt as 

Dtu:=dtu-Ru for R:=dtMM-\ (4.10) 

with M = KS/^, where K, are as dehned in (m and (II.3p . It is easily to be 
known that Df preserves the div^ -free condition, since 

Jdiv^(At’) = J diVj^{MdtiM~^v)) = div{dt{M~^v)) = dtdiv{M~^v) = dt{J divj^v), 

where the equality J div^ v = div(M“^u) can be found in Page 299 of [6]. 

Definition 4.4. Suppose that the forcing functions satisfy 

F^ e L2([o,T];H'i(H))nC°([0,r];H'°(H)), 

F^ e L2([o,r];H'i(H))nC°([0,r];H'°(H)), 

F^eL2([o,T];Fi(S))nC°([0,r];i7^(E)), 

dt{F^ - F^) € L‘^i[0,T]-ioH\n))*), dt{F^ + F^) € L‘^i[0,T];ioH\n))*). 

We also assume that uq G n .^(0) and 9o G n If there exists a pair 

{u,p,9) achieving the initial data uq, 9q and satisfies 

u G L\[0,T]-H^)nC^{[0,T];H^)r] dtu G L'^{[0,T]; H^) n C^{[0,T]; H^) 

9 G T2([o,T];H'3)nC°([0,T];H'2) dt9 G 

QlQ e 

such that they satisfies dni) in the strong sense, we call it a strong solution. 

Then, we have to prove the lower regularity of strong solutions. 

Theorem 4.5. Suppose that the forcing terms and the initial data satisfy the condi¬ 
tion in Definition \4.4\ and that uq, F^(0) satisfy the compatibility condition 

Ho (T^(0) + n^oWiff) = 0, where JpQ = (-9i?7o,-52%, 1), 


(4.13) 
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and Ho is an orthogonal projection onto the tangent space of the surface {x^ = r^o} 
defined by 

Uov = v-{v^o)^o\^o\~^- (4.14) 

Then there exists a strong solution {u,p,6) satisfying ()4.12l) . Moreover, 

+ \\dtU\\\2Hl + \\d'ju\\(^^^Y + HpH^oo/^I + 11^11^2/^2 
+ + II^IIl 2//3 + ||'9t^*|lioo//o + \\dtO\\\2H^ + \\dt^\\{jef)* 

^ -P(II%IIh5/2) (1 + ^iv)) exp {C{1 + J(f{ri))T) (^||uo ||^2 + ||^o ||//2 + ||l^^(0)||^o (4.15) 

+ ll^^(0)ll//0 + 11 (0)11 ^1/2 (J,) + ||F^||i2//l + ||F^||i2//l + ||i^^||^2//3/2(s) 

+ ll'""lli.«./.(E) + + ^■')ll«).). 

where C is a constant independent of r] and (r/) is defined as 

jr(7/):= sup (||r?||^ 9/2 + ||5tT/||^7/2 + ||'9i%||^5/2) • (4.16) 

0<t<T 

The initial pressure, p(0) € H^{Q) is determined by terms uq, 9q, F^{0), F^{0) as a 
weak solution to 

' div^o (V^(,p(0) - F^(0) - 6»oV^o?/3,o) = - div^o(i?(0)uo) G 
< p(0) = iF\0) + B^,uo^o) • ^o\FKo\-^ G (4.17) 

(V^/oP(0) - F^{0)) ■ V = A^qUo • V G i7"^/^(Sfe), 
where y^fl in terms ofr]Q. Also, <9*0(0) satisfies 

dtOiO) = A^„0o + F^{0) G H\n), (4.18) 


and Dtu{0) = dtu{0) — R{0)uo satisfies 

Dtu{0) = A^oUo - V^op(O) + F^{0) + 0063 - 9?(0)uo G ^(0). (4.19) 

Moreover, dtO satisfies 

' dt{dt9)-A^{dte) = dtF^ + G^ in n, 

< V^{dt9)-^ + dt9\M^\ = dtF^ + G^ on S, (4.20) 

_ 5*0 = 0 on Sfe, 


and DfU satisfies 

' dt{Dtu) - A^(Au) + VAdtp) - Dt{9V^y^) = DtF^ + G^ 
div ^(Dtu) = 0 

' SAdtP,Dtu)^ = dtF^ + G^ 

DfU = 0 


in 

in 

on 


n, 

n, 

^ (4-21) 


on Tib, 


in the weak sense of (BAD, where G^ is defined by 

G^ = —{R+dtJK)Aj^u—dtRu+{dtJK+R+R^)V_„/p+divj,/{Oj,/{Ru)—ROj^u+I])Q^j,/u) 

(R~^ denoting the matrix transpose of R), G^ by 

G^ = -dtJKA^9 + div^(-i7V^0 + Va,^0), 


by 


G^ = Bj,/{Ru)M" - {pi - B^u)dtFT + 
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and by 

= -v^e • dt^ - V9,^9 -yK-edt \^\. 

More precisely, (|4.20l) and (14.211) hold in the weak sense of (|4.3I) in that 

+ {dtO, (f) + {dt6\M^\ ,<P)l‘2ho{t) 

= (8,(F’ + + (d,JKF\ 4 >)^ - {dtJKBA 4 >)m« (4,22) 

~ f f ■'S/dtJ^(t>) J 

Jo Jn 

and 

= {dt{F^ - + (5t(6»V^y3),V’)^o - {dtRu + R^tu,^|J)^o 

+ [dtJKF^,^})^0 - {dtJK6e3,'4})j^o - {dtJKdtu,^;)jg,o - (p, 

f [ (dtJKBj^u : : Bg^^V’) J 

^ Jo Jn 

for all (p £ Ip G ^T- 


Proof. Here we will use the Galerkin method, which may be referred to [5]. 

Step 1. The construction of approximate solutions for 6. Since the scalar-valued 
space noH^(H) is separable, we can choose a countable basis Note 

that this basis is time-independent. Now, we need to construct a time-dependent 
basis for F[^ n We define jJ = 4J{t) := K{t)wF According to the Proposition 
12.21 4j{t) G n JF^{t), and {<A^(i)}^i is a basis of n Jif^{t) for each 

t G [0, T]. Moreover, 

dt(tP{t) = dtK{t)w^ = dtKJKw^ = dtKJ(tP{t), (4.24) 

which allows us to express dt(j>' in terms of (jP. For any integer m > 1, we define the 
finite-dimensional space J^Ji{t) := span {(/)^(f),..., (p^{t)} C n J^^{t) and we 

define ^m{t) for orthogonal projection onto J^j^{t). Clearly, 

if 0 G n JF^{t), ^1^6 —)■ 0 as m —oo. 

For each m > 1, we define an approximate solution 

9^ = dJp{t)(fj{t), with : [0,T] —)• M for j = l,...,m, 

where as usual we use the Einstein convention of summation of the repeated index j. 
We want to choose dj* such that 

p)j^o + (d™, + {9"^ 1^1, 4>)h0{T,) = (-^^5 4>) j^o + {F^, 4>) h0(t.) ’ 

with the initial data d™'(0) = JPf^9Q G for each p G J^jp{t). And (I4.25P is 

equivalent to the system of ODEs for 


“j 



+ d- 




JifO 


+ 



for j,k = l,...,m. The m xm matrix with j,k entry is invertible, 

the coefficients of the linear system (j4.26p are T]), and the forcing terms are 
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cO([o,r]), so the usual well-posedness of ODEs guarantees that the existence of a 
unique solution d™ G (^^([OjT]) to (I4.26P that satishes the initial data. This provides 
the desired solution, 9^, to ()4.25p . Since F^, satisfy (14.lip , equation (14.261) may 
be differentiated in time to see that d™ G T]), which means d™ is twice 

differentiable almost everywhere in [0, T]. 

Step 2. The energy estimates for d™. Since d™'(t) G we take 0 = d™' as a 

test function in (j4.25p . using the Poincare-type inequalities in Lemma A. 14 of [6] and 
usual trace theory, we have 

< i\\F\^, + ^ |d™|20,j. 

Then, applying Cauchy’s inequality, we may derive that 

with Co{r]) := I + supo<t<T 11 <9* <^-^11 Using the Lemma 2.9 in m, we may have 
||d™(0)||^o < Pi\\vo\\H^/ 2 )\\e^io)\\HO < Pi\\r]o\\H^/^)\\e^io)\\m 
= -P(lldollH5/2)||<^o^do||//2 < P{\\'no\\H^/2)\\9o\\H2- 


Now, we can utilize Gronwall’s lemma to deduce energy estimates for d”*: 

sup \\e-\\%o + \\e-\\U 

0<t<T T 

< P(||r?o||H5/2)exp(Co(7?)r)(||do||^. + \\F^\\%o + 


Step 3. Estimates for dt9^{0). If d G H‘^{Q) (b 4> £ the integration by 

parts reveals that 


(d, (p)— j —A^6(j)J F j J2/9 ■ ^)4> — i—Aj2/6, + (V^d • c/U, (P)ho{t.) 

(4.29) 

Evaluating (I4.25P at t = 0 and employing (I4.29p . we have that 

{dt9^{0), cP)^o = (A^o^™(0) + i""(0), <^)^o , (4.30) 

for all (j) G 

By virtue of ()4.24p . we have that 

dtO^ - dtK{t)j{t)e^{t) = G jCit), (4.31) 

so that (p = 5td™’(0) — (9tAr(0) J(0)d”’’(0) G J^(0) is a choice for the test function in 
(I4.30p . So using this test function in (I4.30p , we have 

||9,d-(0)||^o < ||a,iL(0)J(0)d-(0)||^o||9id™(0)||^o 

+ ||5td™(0)-atA:(0)J(0)d™(0)||^o||A^„d-(0) + F3(0)||^o. 


Then after using (14.271) and Cauchy’s inequality for the right-hand side of (14.321) . we 
have the bound 


\\D,rm\F S. GW (l|9olli> + I|F’(0)|I5^„) (4.33) 

with Ci(ri) = P{||))(]||„vs) (1 + l|8iA'(0).4(0)||J,« + UMifci)- 

Step 4. Energy estimates for dt9^. Now, suppose that (p{t) = FJ^{t)<pP for cT G 
C°4 Qo,T]), j = 1,..., m; it is proved as in (I4.3ip . that dt(p — dtK{t)J{t)(p G =^^(t) 
as well. Then in (j4.25p . using this (p, and temporally differentiating the result 
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equation, and then subtracting from the result equation (I4.25P with test function 
dtcf) — dtK{t)J{t)<j), we find that 


(020-, + (5*0- 1^1 , <^)^0(S) 

= (5*(f3^f5),</.)(^i). + (F3,(5*irj + 5*JK)</.)^0 + (F^5*KJ</.)^0(^^ 

- (5*0-, (5*i^J + dtJK)(j))^o - (0-, dtKJ4>)^, - (0-, 


- [ (5* Ji^V^0-• V^(/) + Va,^0-• + V^0" 

Jn 




J. 


According to (14.311) and the fact that d—(t) is twice differentiable almost everwhere 
as we have pointed in the first step, we use (j) = 5*0— — dtKJO'^ as a test function 
in (|4.34p . Utilizing Cauchy’s inequality, trace theory and the Remark 2.3 in [6], we 
have that 


dt - (5*0-,5*a:j0-)^o) + 

< Coiri) (^||0”^||^o - (5*0-,5*A'J0-)^o) +C2(r?)||5*0-||^* (4-35) 

+ C’ (llF'll^o + + C||5*(f3 + 


for C 2 {r]) is defined as 

C2{v) ■ = sup [l + ||5*(5*A:j)||ioc + ||5 *A:J||^i + ||5*i2/||ioc 

0<t<T 

+ (1 + ||32/|||oo)(l + 115 * 7 x 11100 )] (1 + ||5*XJ||^i). 

Then according to Cauchy’s inequality and Gronwall’s lemma, (j4.35p implies that 

sup (||a*0-||^o + ||9*0-||^i 
0<t<T 


exp(C„(„)T)(|| 8 , 9 "'( 0 )|| 3 ^. + C,(,)||()”( 0 )|| 3 ^. + WF^Mf 

+ I|r^|li2;,v2 + + r^)ll^^^i).) 

+ C 2 {r])( sup ||0-||^o+ / exp(C'o(? 7 )(r - s)) 

\ 0 <t<T Jo 


(4.36) 


(«) 


ds 


Now, the energy estimates for dtO^ is deduced by combining (I4.36P with the estimates 
(ICTll . (ICTIl and (031) . 


sup ||5*0-||^o + ||5*0-||;^i 

0<t<T 


< 


(Ci(r?) + C 2 (r?))exp(Co(r?)r) (||0-(O)||^o + ||F'(0 )||^o) 


(4.37) 


+ exp(Co(r7)r) |C2(r?) [\\F^\?^o + ||T^|| )+mF^+F^)\\l^r^ 


i5||2 


?5m|2 
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Step 5. Improved estimates for 0'". Using the 4> = dtO^ — dtKJO^ G as a 

test function in (I4.25p . we can improve the energy estimates for 0™. 

= {dte^, dtKje^)^o + ( 0 "*, dtKje^)^^ + {f^ dtO^ - dtKje^)^, ( 4 . 33 ) 

+ {F\dte^ - + 1^ (v^e^ • ?j . 


Since we have already controlled ||^™'||^i and ||5t0™'||^i, integrating (I4.38P in time 
implies that 


sup \\o^rj^, + \\dt9^r^, 

0<t<T ^ 


< 


P{\\r]o\\H^/2) {C,{r^) + C2{v)) exp(Co(ry)r) (||0o||^o + ||F'(0 )||^o) 


+ P{\\Vo\\H^/2)exp{Coiv)T) C2{v) [WFYj^o + ||i"1IW/2 




(4.39) 


Step 6. Uniform bounds for (|4.37p and (I4.39p . Now, we seek to estimate the 
constants Ci{r]), i = 0,1,2 in terms of the quantity JF^rj). A direct computation 
combining with the Lemma A. 10 in [6] reveal that 


Coiv) + Ci{rj) + C 2 {ri) < ^(1 + (4.40) 

For a constant C independent of rj. 

Step 7. Passing to the limit. According to the energy estimates p4.37p and p4.39p 
and Lemma [2.11 we know that the sequence {9”^} is uniformly bounded in L°°H^ and 
{dt9^} is uniformly bounded in L°°F[^r\L^H^. Then, up to extracting a subsequence, 
we know that 


^ 0 weakly- * inL°°F7^, dt9^ 9f0inL°°F7°, ^ weakly inL^Fl^, 


as m —>• oo. By lower semicontinuity, the energy estimates reveal that 

+ \\^t9\\\^H° + \\^t9\\\2Hi 

is bounded from above by the right-hand side of p4.15p . 

According these convergence results, we can integrate p4.34p termporally from 0 to 
T and let m ^ oo to deduce that d^9'^ 8^9 weakly in with an action of 8^9 

on an element (j) G defined by replacing 0™ with 9 everywhere in (14.341) . From 
passing to the limit in (I4.34p . it is straightforward to show that is bounded 

from above by the right-hand side of (I4.15p . This bound shows that 8t9 G C^L^. 

Step 8. In the limit, (I4.25P implies that for almost every t, 

{8t9,(j))j^o+i9,(j))j^i+i9\J^\ = {F^,(j))^o + {P^,4 >)for every cj) G . 

(4.41) 

For almost every t G [0,r], 9{t) is the unique weak solution to the elliptic problem 
P3.ip in the sense of (I3.2p . with F^ replaced by F^{t) — 8t9{t) and F^ replaced by 
F^{t). Since F^{t) — 8t9{t) G F7°(ll) and F^{t) G 77^/^ (S), Lemma lOl shows that 
this elliptic problem admits a unique strong solution, which must coincide with the 
weak solution. Then applying Proposition 13.51 we have the bound 

l|9(i)llL < C(.,„) (||8,«(t)ll5f.-> + \\FHmirr-, + l|F"(*)ll?,,-,„(,;|) (4.42) 
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when r = 2,3. When r = 2, we take the superemum of (I4.42p over t G [0, T], and when 
r = 3, we integrate over [0, T]; the resulting inequalities imply that 6 G 
with estimates as in (I4.15p . 

Then for the linear Navier-Stokes equations, the process is exactly the same as 
[6]. Then we know that {u,p, 9) is a strong solution of (jl.Sp with the estimates as in 

dnn]). 

Step 9. The weak solution satisfied by dt6 and Dtu. We may integrate (I4.34p in 
time from 0 to T and pass the limit m —)• oo. For any (j) G we have dtKJcj) G 
so that we may subsititute dtKJcp for cj) in (14.411) : this yields 


{d^9, + {dtO, 4>)^i + {dtO \^\, 4 ‘)l^ho{j:) 

= {dt{F^ + + {dtJKF\4>)^, - {dtJKdte,cf)^o 

~ f f ji/9 ■ S/si/4> FJ 

Jo Jn 


(4.43) 


for all 4> G This is exactly the (I4.22p . To justify that (14.221) implies (I4.20p . we 
may integrate by parts for the equality 


f [ {dtJKS/^O ■ S/^(f) + gts/9 ■ 'Vsi/cj) + ^s/9 ■ Vg, 

Jo Jn 

l>T f 


J 


(4.44) 


10 Jn 


= {div^{-R'V^u + Vg^^u),^)jg>o - {V^u-dt^ + ^dts/u ■ ^A) • 


We then may deduce from p4.22l) that dt9 \s a weak solution of (14.201) in the sense of 
(14.31) with dt9{R) G Then we may appeal to the computation in [6] to deduce 

that p(0) satisfies the equation (I4.17P and Dtu is a weak solution of (I4.2ip in the 
sense of (|4.3p with Dtu{0) £ &'{0). □ 


4.3. Higher regularity. In order to state our higher regularity results for (|1.8I) . we 
need to construct the initial data and compatible conditions. First, we define the 
vector or scalar fields in 12 and on S by 

<t^^{G^,v,q) = A^u-+ - Rv, 

(£°2(g3,0) = A^0 + G^ 

£^(u, q) = -{R + dtJK)A^v - dtRv + {dtJK + R + 

+ div^(D^(i2u) - RBj^v + Og^^v), (4.45) 

^3(0) = -dtJKA^e + div^(-i?V^0 + Va,^0), 
e^(u, q) = B^(i?u)^ - {ql - 0^v)dt^ + 

^5(0) = _V^0 • dt^ - Vg,^e -^-Qdt \FK\ , 

and we define functions in 12, on S and on by 
f^(G\u) = div^(G^ - Rv), 
fiG\v) = (G^ + O^v^) ■ 
f{G\v) = iG^ + A^v)-iy. 


(4.46) 
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We write + ^V^ys, F^’^ = F^, F^’O = F^ and F^-O = F^. When F\ 

F3, F^, F^ u, p, and 6 are regularly enough, we can recursively define 

Fhi - DtF^’^-^ - d{-\eV^y3) + Di-\eV^y3) + 9^-^) 

= DiF^ - [d{-\9V^ys) - Di-\9V^y3)) + ^ 

e=o (4.47) 

F^J := + <B^{d{~^9) = d{F^ + J2df<B^idi~^~^9), 

e=o 


in n and 




F"^’^ := dtF^’^ ^ + e:^(Fj ^u,d{ ^p) = dlF'^+ ^ ^u,d{ ^ ^p), 


£=0 


i-i 


f 5J — a^irs-i-i + (B^{d{-^9) = d{F^ + ^ df(B^{d{-‘^-^9) 


e=o 


on S, for j = 1,..., A^. 

Now, we define the sums of norms with F^, F^, F^ and F^. 


Af-l 

5(F^F^F^F5) := ^ 

j=0 

+ \\d^ F^\\L2(oH^{n))* + \\d^ -^^IIl2(ph1(o)) 

N 

+ E 

j=0 
N-l 

+ E 

i=o 

Af-l 

+ \\d{F^ 


\]^2fj2N-2j-l + ||5jF^||^2j:^2JV-2i-l 

'^'|Il2(c 

\\L^H^N-2j-l/2 + \\d{F^\\ ]^2jj2N-2j-l/2 
J^cofj2N-2j-2 + \\dlF'^\\j^oofj2N-2j-2 


+ E 

i=o 

Af-l 


(4.48) 


(4.49) 


||l,oojy2JV-2i-3/2 


lcofj2N-2j-3/2 I , 


:= E (l|9/i"'(0)||H2iv-2.-2 + ||a^F3(0)||^2.-2,- 

j=0 

N-l 

+ E 

i=o 


jj2N-2j-3/2 + ||l9t F^(0)||^2iV-2j-3/2 


For simplicity, we will write 3^ for 3^(F^, F^, F^, F®) and Jo for ^q{F^,F^,F^,F^) 
throughout the rest of this paper. From the Lemma A.4 and Lemma 2.4 of [6], we 
know that if J < oo, then 

diF^ G c'°([o,r];F2^-2j-2(Li)), d{F^ G c^{[o,Ty, 

d{F^ G C'°([0,r];F2^-2j-3/2(s)), and d^F^ G C°([0, T]; F2^-2J-^/2(S)) 
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for j = 0,..., N — 1. For ij, we define 

Af+l 

'^{v) •= \\d{rj\\'^j^2j^2N-2j+5/2, 

i =2 

N 

<^iv) ■= WvWl cx>^ 2 iV+l/ 2 (s) + Y.\\9 {v \\1 aoH'2N-2j+3/2(^Y.y 

i=i 

^{r]) := D{r]) + ^{rj), 

N 

^ oiv ) •= ||%||^ 2 iV+l/ 2 ( 5 ^) + ||^^( 0 )||^ 2 JV- 2 j+ 3 / 2 (-y;V 

i=i 

These following lemmas are similar to Lemma 4.5, 4.6, 4.7 in [6] as well as the idea 
of proof, so we omit these details here. 


(4.50) 


Lemma 4.6. If k = 0,..., 2N — 1 and v, 0 are sufficiently regular, then 

\\dtv - DtvWl^jjk < P{Ii{ri))\\v\\l 2 Hk^ 

Ii 9 i(ev^y 3 ) - A( 0 v^i/ 3 )iii 2 ^. < 
and if k = 0,..., 2N — 2, then 

\\dtv - DtvWl^^k < P{K{'n))\\v\\l„^jjk, 

mev^ys) - A(0v^y3)iiiooH. < PimmwlooHk- 
If m = 1,..., N — 1, j = 1,... ,m, and v, 0 are sufficiently regular, then 

\v\\i^ 2 H'im- 2 j +3 + \\d^v\\^^^ 27 n- 2 j +2 j , 


(4.51) 

(4.52) 

(4.53) 

(4.54) 


iv - Dlv\\l2H2m-22+3 < P{Ii{v)) (|| 5 , 


1=0 


i-1 


\\div - Div\\l„,jj2m-2j+2 < P{Ii{v)) Y \\dtV\\YH^m-2j+2, 


(4.55) 

(4.56) 


e=o 


\\d{{&Vji/y3)\\'j^2H2m-2j+2 fS P{K{r]))Y^ (^\\df&\\‘j^2jj2m-2j+3 + ||5f0||^ooJ|^2m-2j+2 

£=0 

(4.57) 

Wdiiev^ys) - Di{QV^y5)\\Y^2m-2j+3 < P{m) E 11 ® 11 XjOafj2m-2j + 2 , (4.58) 

i=0 

and 

\\a,DTv - + \\i)}DTv - arMlx-r)- 


V £=0 


>H2 


(4.59) 


Also, if j = 0,... ,N and v is sufficiently regular, then 

llc^KO) - O^( 0 )||^„-„ < P(€„{ri)) Y, l| 8 fi>{ 0 )||?,„-=,, 


£=0 


(4.60) 
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and if j = 0,..., N — 1 and 0 is sufficiently regular, then 

\\d{{e{0)V^,y3,o) - Dl{e{0)V^,y3^o)fH2N-2,-2 < P{Mv)) E 115^0(0)11^2.-2.-2. 

£=0 

(4.61) 

Here all of the P{-) are polynomial, allowed to be changed from line to line. 


Lemma 4.7. For m = 1,... , N — 1 and j = 1,... ,m, the following estimates hold 
whenever the right-hand sides are finite: 


I II 2 II Tp3,j ||2 

\r \\ 1^2fj2m-2j+l ~r Hr 


l2^2m-2j+l 


_1_ II ||2 I II ||2 

— 2.7 + 1 “T ||.e \\^2^2m.—2j+3/2 “T ||.e ||r 


j^2jj2m—2j+3/2 


i-1 


< 


P{m)h + E(ii dt'^\\‘L2 H2m-2e+3 + II 5f0 11^2/^2771-2^+3^ 

^ £=0 

+ E (l|5t^llioo//2m-2£ + 2 + ||5f0||^oo//2m-2£ + 2 + 11 5fp 11 ^2//2m- 

£=0 

+ l|5fp||EH2m-2£+l^ j , 


(4.62) 


p^d 


2 

/_oo//2m-2. 


+ IIF^J 


2 

]^oo fj2m — 2j 


+ ||F^’^' 


||2 

11 Jj2m—2j-\-l/2 


+ \\F^’^ 


||2 

11 £jOo Jj2m—22 


< P{R{g))(^ +E (\\dtU\\‘j^ooJj27n-2e+2 + ||5f0||E//2m-2f+2 

^ £=0 


(4.63) 


+ l|5fp||E//2m-2f+l^ J , 

_ F^n\\hi,HHn)r + l|5i(i^^’”^ + i^^’”^)lli2(„//i(^))* 
< P{A{r]))^ + ||5™t(|||2//2 + ||57*6'|||2//2 + ||57 *p|||2//i 

m—1 

+ E (l|5t^llL°oJ?2 + ||5fu|||2//2 + ||5f0||E//2 + ||5f0|||2//2 

£=0 


+ l|5fp||ioo//l + ||5fp||^2//2^^ . 

Similarly, for j = 1,..., N — 1, 

11 .^^’■^'(0)11 ^2.-2.-2 + ||1^^5(0)||^21V-2.-2 + 11 (0) 11 ^2JV-2. -3/2 + 11 (0) 11 ^2JV-2j-3/2 

< P((2:o(?/))( 5^0 + ||5/0(0)II//2JV-2/ + E (l|5/^f(0)||//2.-2£ m 

e=Q ^' 

+ l|5f0(O)||//2JV-2£ + ||9fp(0)||//2JV-2£-l)^. 


Here all of the P{-) are polynomial allowed to he changed from line to line. 

Lemma 4.8. Suppose that v, q, G^, are evaluated at t = 0 and are sufficiently 
regular for the right-hand sides of the following estimates to make sense. If j = 
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0,..., — 1, then 

\\<B^^{G^,v,q)\\]^2N-2j-2 

, o n 1 o X (4.00) 

<^’(e^o(r/))(|| v\\jj2N-2j + ||(?||^2JV-2j-l + I|G'|| ‘]j2N-2j-2) , 

0)||^2JV-2j-2 < P{^o{r])) (||0||^2Ar_2j + 11 11 ^2JV-2j-2 ) ■ (4-67) 

If j = 0,..., N — 2, then 

\\f{G^,v)\\H2N-2i-3 + ||f^(G^,i;)||^2JV-2i-3/2 + 11 (G^ , II) 11 ^2JV-2i-5/2 

^ Pi^oiv)) (||G^||^2Ar_2j-2 + ||G'^||^2JV-2i-3/2 + 11II11 ^2JV-2i ) • 

For j = N — 1, i/div^(o) ii(0) = 0 in Q, then 

llf^(G^,i;)||^i/2 + ||f^(G^,i;)||^_i/2 < P{^o{r])) (||G^||^2 + ||G'^||^i/2 + ||ii||^2) . 

(4.69) 

Here all of the P{-) are polynomial allowed to be changed from line to line. 


Now we can construct the initial data and compatible conditions. We assume that 
uq e 6 q S , r]Q S Then we will iteratively construct the 

initial data d{u{ 0 ), d{ 9 { 0 ) for j = 1 ,..., and &ip{ 0 ) for j = 1 ,..., N" — 1 . First, 
we denote FbO(O) = e = F^{ 0 ) G F^’^iO) = F^{ 0 ) G 

pp 2 N- 3 / 2 ^ F 5 ’ 0 ( 0 ) = F^{ 0 ) G 772 ^- 3/2 and D^uiO) = uq G 9 ?^( 0 ) = Oq € 

. Suppose now that we have constructed G 7727V-2£-2^ p 3 ,e ^ j^2N-2i-2^ 

pA,e g p2N-2l-3/2^ p 5 / g ^2V-2£-3/2^ G 9 ^ 0 ( 0 ) G for 

0 < 7 < j < A^ — 2 ; we will construct d{p{ 0 ) G 772^-27-1 g^g gg G 

jj 2 N- 2 j- 2 ^ a^+i 0 (O) G 772^-27-2^ FiJ+i( 0 ) G 772^-27-4^ F 3 J+ 1 ( 0 ) G 772^-27-4^ 
^4,7+1 (0) g jp2N-2j-7/2 gg^ 7 ^^’ 7 + 1 ( 0 ) G 772'^^-27-7/2 g^g follows. 

By virtue of estimate, we know that 

/I = f\F^’\ 0 ),D{ui 0 )) G 772^-27-3^ 

/2 = f 2 (F‘^d(o), 77 ^ii( 0 )) G 772^^-27-3/2^ (4 

/3 = f (f 1 - 7 (o), 77 ^u( 0 )) g 772^-27-5/2 

This allows us to define c^p(O) as the solution to ( 13 . 181 ) . The choice of /^, / 2 , / 3 j 
implies that dfpfS) G 772^“27-i^ according to the Proposition 2.15 of [T 4 j. Now the 
estimates (I 4 . 65 jl . (Id.GOjl and (I 4 . 66 jl allows us to define 

7 ?^'+'u( 0 ) := eioi + 5 ^( 0 (O)V^oy 3 ,o),^^ii(O), 5 ^p(O)) G 

d{^^ 0 { 0 ) := (^F 3 d(o), 5 ^ 0 ( 0 )) G 772^-27-2^ 

77iJ+i(0) := DiF^’^O) - di{e{0)V^,y3,o) + Di{e{0)V^,y3,o) 

+ ( 77 ^u( 0 ),a^p( 0 )) G 772^-27-4^ 

^3,7+1 (0) atF 3 j( 0 ) + (^d{e{ 0 )^ G 772^-27-4^ 

F^- 7 + 1 (o) := (^D{u{ 0 ),d{p{ 0 )^ G 772^-27-7/2^ 

^ 5 , 7 + 1 ( 0 ) := 5 t 77 ' 5 ’^( 0 ) + (5^0(0)) G 772^-27-7/2^ 
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Then, from the above analysis, we can iteratively construct all of the desired data 
except for D^u{0), d^~^p{0) and d[^6{0). 

By construction, the initial data Dlu{0), dlp{0) and d^9{0) are determined in terms 
of uq, 00 as well as dfF^{0), dfF^{0), dfF‘^{0) and dfF^{0) for £ = 0,. .. , N — 1. In 
order to use these in Theorem 14.51 and to construct D^u{0), d^~^p{0) and d[^9{0), 
we must enforce compatibility conditions for j = 0,...,iV — 1. We say that the j-th 
compatibility condition is satisfied if 

[ Diu{0) G 

I Bo + B^,d{u{0)A>) = 0. 

The construction of Dlu{0) and dlp{0) ensures that Dlu{0) G F['^{f}) and div^g(T)ju(0)) 

0. 

In the following, we define d[^0{O) G Fl^, d^~^p{0) G Fl^ and D^u{0) G F[^. First, 
we can define 

d^9{0) = (a°^{F^’^-\0),d^-^9{0)) G 

employing ()4.67p for the inclusion in Then using the same analysis in [6], the 
data df~^p{0) G F[^ can be defined as a weak solution to (I3.18p . Then we define 

Dfu{0) = [f^^^-\ 0) + df-\9{0)V^,y3,o),D^-\{0),df-^piO)) G 

employing (I4.66P and p4.6ip for the inclusion in H^. And D^u{0) G ^(0) is guar¬ 
anteed by the construction of dj^~^p(0). Combining the inclusions above with the 
bounds (I4.65p . (I4.68h , (I4.66h and MTTh implies that 

N N-1 N 

||Tlj n(0)||^2JV-2j + ^ \\d{p{0)\\‘j^2N-2j-l + ^ Il'St ^(0)||^2JV-2i 

j=0 j=0 j=0 (4-72) 

< Pi^oir])) {\\uo\\ H2N + ||^o||^2JV + 5o) • 

Before stating the result on higher regularity for solutions to pi.8p , we define some 
quantities: 


D{u,p,9) 


(£(u,p,6») 

^u,p,9) 


N 

^ (ll^^lli2j^2JV-2j+l + \\d{9\\‘j^2H2N-2j + l'^ + (SV)* 

j=0 

N-1 

1=0 


(4.73) 


N 

3=0 




:= ^{u,p,9) + (B{u,p,0). 


N-l 

3=0 


Theorem 4.9. Suppose that uq G 9o G po G and 

d <oo. Let d{u{0) G di9{0) G and d{p{0) G 

for j = 1,N — 1 along with D^u{0) G (0) and d[^9{0) G F[^, all be determined 
in terms of uq, 9o and d{F^{0), c^F^(O), d{F^{0), d{F^{0) for j = 0,... ,N — 1. 
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There exists a universal eonstant Tq > 0 such that ifO<T <Tq, then there exists 
a unique strong solution {u,p,6) on [0, T] such that 

4u G CO n L2 ([0,r];/72Af-2i+i(O)) forj = 0,...,N, 

d{p G CO {[0,T]-H^^-^^-\n)) n l2 {[0,Ty,H^^-^^{n)) forj = 0,...,N-1, 

die G ([ 0 ,r];F 2 ^- 2 ^(O)) n {[0,T]-,H‘^^-^^+\n)) forj = 0,...,N, 

G {^tT, and df^+^e G . 

The pair {Diu, dip, di9) satisfies 

dt{Diu) - A^{Diu) + V^{dip) - di{OVj^y^) = inQ, 

diy j^{Diu) = 0 mfi, 

dtidie) - A^idie) = inn, 

S^idip, Diu)FT = F^'^ onT, 

V^{die)-F^ + die\F^\=F^^^ onT, 

Diu = 0, die = 0 onTib, 

in the strong sense with initial data ^Dj u(O), 5/p(0), (9^ 0(0)^ for j = 0,... , N — 1, 

and in the weak sense with initial data D^u{0) G ^F{0) and d^e{0) G H^. Here the 
forcing terms F^’^, F^d ^ and p^d are as defined by (I4.47P and (I4.48p . Moreover, 

the solution satisfies the estimate 

^{u,p,e) < p(£o(??),-<^(?/))exp (rp(£(77))) {\\uo\\jj 2 N + \\eo\\H 2 N +do + d), ( 4 . 75 ) 

where the constant C > 0, is independent of rj. 


(4.74) 


Proof First, notice that P{-,-) and P(-) throughout this proof is allowed to change 
from line to line. Theorem 14.51 guarantees the existence of {u,p,e) satisfying the 
inclusions (I4.12p . The {Diu, dip,die) are solutions of (I4.74p in the strong sense when 
j = 0 and in the weak sense when j = 1. Finally, the estimate ()4.15l) holds. 

For an integer m > 0, let Pm denote the proposition asserting the following three 
statements. First, {Diu,dip,die) are solutions of (I4.74p in the strong sense for j = 
0,..., m and in the weak sense when j = m + 1. Second, 

diu G L°°i 72 "*- 2 j +2 p J2jj2m-2j+?.^ qJq ^ ^oo^2m-2j+2 ^ j2^2m-2j+‘S 

for j = 0,1,..., m + 1, d^^'^u G (^r)*j df^^‘^9 G and 

dip G n piH‘im-2j+2 

for j = 0,1,..., m. Third, the estimate 

m+l 

{\Wt'a\\\oofj2m-2j+2 + \\diu\\\2H2vn.-2j+3 + 11 ^ 11 ^oo j:^2m-2i+2 + 11 6* 11 ^2 j:^2m-2i+3 ^ 

J=0 

+ 4" {\WtP\\L‘^ H'^-'n-23 + l + \\dip\\\2H2m-2j+'^ 

j=0 

< P(£o(h),-^(?/))exp {TP{(B{p))) (||uo||^ 2 iv + \\eo\\jj 2 N +50+5) 

holds. 
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We will use a finite induction method to prove that Pm holds. Theorem 14.51 implies 
that Pq holds. Then in the rest of this proof, we will divide the proof into two steps. 
Step 1. Proving the first assertion. Suppose that Pm holds for m = 0,..., Ai — 2. 
From (j4.62l) - (l4.64l) of Lemma W7A we have that 

( m 

d + ^ {\\dtV\\‘j^2H3 + ||9f0|||2j^3^ (4.77) 

e=o 

m 

+ ^ + \\dfq\\‘^j^2H2 + \\dfq\\‘^ioo 

e=o 


y m 


(4.78) 


+ l| 9 ,(FS.’"+i(e) - F^’"+i(e))lli=, 

<p(ji(.,))(5 + ||er'f|l + liar'eili^H- + liarVll 

m 

+ ^ m + ll^t'^^lli2_H-2 + ||9f0|||oojy2 + 1119^01112^^2 

£=0 

+ + l|5t9llL2j|^2^^ . 


(4.79) 


Now we will use the iteration method. We let vP be the extension of the initial 
data (9^u(0), j = 1,..., A^, given by Lemma A.5 in [6], which may also give 9^, the 
extension of the initial data (9^0(0), j = 1,..., A^, and similarly let be the extension 
of d(p{0), j = 1,..., N — 1, given by Lemma A.6 in [6]. By (j4.72l) and the estimates 
given in the Lemma A.5 and Lemma A.6 in [6], we have 


N 


d 0||2 


{\\diu^\\l2H2N-2j+l + \\diu^\\l^^2N-2j + ||5^0°|| \2jj2N-2j+l + HO' I \L^H2N-2j 

j=0 

N-l 


+ E 

j=0 


. 0||2 


., 0||2 


N 


'tl' \\L‘^H‘^N-2j + ||C^P' ||j;^cx3J|^2JV-22-1 

N-1 


(4.80) 


N 


^ ^ ||99j W(0)||^2iV-2j + ^ ||^p(0)||^2]V-2j- 1 + ^ ||c^6*(0)||^2iV-2j 

i=o i=o j=o 


Pi < Bo { v )) (lluoll H 2 N + ll^0||^2JV + T?o) 








BENARD CONVECTION 


33 


According to (I4.77p - (l4.80p . we may derive that ,p^), {u^, 

and F^’'^~^^{9^) satisfy (|4.11l) . Also the compatibility condition (I4.13P with F'^ re¬ 
placed by and uq replaced by holds by (I4.7ip since and 

p^ achieve the initial data. Then we can apply Theorem 14.51 to find a pair q^,Q^) 
satisfying the conclusions of the theorem. For simplicity, we abbreviate (USD as 
C{v,q,e) = ¥ = {F\F^,F^,F^). Then 


C{v\q\Q^) = := ^^’”^+^( 0 °)), 

vHO) = Dr^^iO), qHO) = 9r+V(0), 0^0) = 5r+'0(O). 

If we denote the left-hand side of (j4.15p as 0), then we may combine (j4.15p . 

(|4.65p . (|4.77p . (|4.79l) and (|4.80p to derive that 

nv\q\e^) < exp (F(e:(0))T) (||uo||^2iv + ||0o|lH2iv + 5o + 5) • 


Now, suppose that {v'^,q^,Q^) is given and satisfies ISPu”, g”, 0”) < oo, we define 
{vF,p^,6^) which satisfies the ODEs 


O/u” (0) = u” (0) for j = 0,..., m. 


(4.81) 


dip^{0) = g”(0) for j = 0,... ,m, 


(4.82) 


r = 0”, 

{ (4.83) 

D c^0"'(O) = 0”(O) for j = 0,..., m. 

From the wellposedness theory of linear ODEs, we know that these ODEs have unique 
solutions. If we define .^(u, g, 0) by 

m 

Si{v,q,B) : = \\dr*'v\\l.„, + l|ar‘ 9 llI>H. + 11 ^*‘ 0 |II>H> + E (ll«l'’h“/f> 

£=0 

+ + \\dt^\\\2fj3 + ||l9f0|||ooj:^2 + 11 g 11 ^2 j:^2 + 11 g 11 |<x) ^ , 

then the solutions of ()4.81l) - ()4.83p satisfy the estimate 

( m 

'F\\aiu(o)f„, + \\4pm?H^ 

(4.84) 

+ ||a^ 0 (O)||^ 3 +r®(u-,g-, 0 -)') <oo, 


where P{T) is a polynomial in T. 

Applying Theorem 14.51 iteratively, we can obtain sequences {(u”, g”, 0”)}^]^ and 
{u"',p”, 0”}^]^ satisfying (|4.81l) - (l4.83p and 

C{v^, q^, Q^) = ¥^+\u^-\p^-\e^-^), 

v^{0) = D^+^iO), g"(0) = Or+Vo), 0*^(0) = 0(0). 


Then 

C{v^+^ - g’"+^ - q'^, 0"+^ - 0") = ¥^+^{u'^ - u'^-\p'^ - p'^-\ 0" - 0’""^), 

u^+HO) - u”(0) = 0, g”+^(0)-g”(0) = 0, 0”+^(O) - 0’^(O) = 0. 
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Since the terms involving F^, F^, F^ and F^ are canceled in — 

pn-i^ Qn _ -vve can use (|4.77l) and (I4.79P to derive that 

||h [U -U ,P -P )\\L2Hi + \\^ [U -U )\\l2hi 

, II 77'4,m+ln „,n—l n—1\||2 ,\\Tp5,m+l(Qn Qn—1\\\2 

+ ||h ’ [U -U ,P -p 111^2^3/2 + ||h [ti -ti 111^2^3/2 

+ ||8,(F‘’”+‘(a” - -P”-') - -p"-‘))lli>,,„i,n|,. 

+ ||8,(F’’”+‘(9” - 9”->) - F^”+‘(«” - «"-‘))lli.(.«.(n))- 

< P(jS(t,)- ii"-', p" - p"-', fl" - 9"-‘). 

Since, for each n, 0”") achieves the same initial data, similar to the ODEs 

()4.8ip - ()4.83l) . we have that 

A{u^-u^-\p^-p'^-\d^-e^-^)<P{Si{p))TP{T)^{v^-v^-\q'^-q'^-\Q^-Q^-^). 

(4.86) 

The above two estimates with (|4.15p imply that 

- v^, - q^, 0’"+^ - 0”) 

<P(C:o(r7),i^(??))exp(P(e(r?))r) 

X rP(r)iB(u” - v^-\q^ - 0” - 0”"^), 


(4.87) 


which implies that there exists a universal Tq > 0 such that if T < Tq, then the 
sequence {(f", 0”)}))T^ converges to {v,q,Q) in the norm y^iB(-, •), which reveals 

that {{u'^,p^, 0"')}^=i converges to {u,p,6) in the norm •)• 

By passing to the limit in (|4.8ip - (l4.83l) . we have that v = q = d'^'^^p and 

0 = +^0. Then, passing to the limit in ()4.85p . we have that 

c{D^+\,dr^^p,dr+^e) = ¥^+\u,p,e). 

Then Theorem 14.51 with the assumption of Pm; which provides that 

d^^^p, (9™^^0) are solutions of ()4.74l) in the strong sense for j = 0,... , m, enables us 

to deduce the first assertion of Pm+i- 

Theorem 14.51 together with the estimates ()4.62p , ()4.79p and (I4.76P , gives us that 
<P(eo(r/),J^(h))exp (P(C:(77))r) (lluoll jj2N + ll^ollpaN (4.88) 

+ 5o + 5+l|8r‘«ll + ii9rviii.„. + ii8r‘«iii>H>)' 

0n the other hand, the estimate (I4.56P implies that 

ii9r+'«ii L2//2 + ii8rviii.„. + iiar‘«iiw 

< T {ii«r+'«lii«p« + li»r‘p|li»p. + liaT'«ili»p4 

< T (iiar+N - or+ llor'«lli~p. + liar Vlli»«. + lisr 

m 

< T(p{m) E + nDT^^u,dT^^p,dT^H)) (4.89) 

1=0 

< T(^P{(Bo{p),^{r]))exp {P{<E{r]))T) (|| Uo\\jj2N + ||^0||^2JV + T?0 + 5") 
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where in the last inequality, we have used (14.761) again. Combining the above two 
estimates, we may further restrict the size of universal Tq > 0 such that if T < Tq, 
then 


< P{<£o{r]),R{r]))exp {P{(B{ri))T) (||uo||^2iv + \\Oo\\h 2 n +5 o + 5')- 


(4.90) 


Step 2. Proving the second and third assertions. In the following, the second and 
third assertions will be derived simultaneously. The estimate of (14.901) with Lemma 
14.61 and estimate (I4.76j) imply that 

/m+2 \ 

< PiMi)) E \\dtU‘\\‘l,2jj2m-2t+3 + \\dfu\\‘^j^aoH2m-2e+2 j 

\^=0 / 

+ P{^o{ll)^ ^iv)) {P{^{'n))P) (||'^^o||^2iV + ||0o||^2JV + T?0 + 5^) (4.91) 

< P{ii{r]))P{^o{rj),R{r])) exp {p{(B{r]))T) {\\uo\\h 2 n + \\Oo\\h 2 n + So + 5^) 

+ P{^o{v)j ^iv)) {P{^{'n))P) (ll^o||^ 2 iv + II^oIIh^jv + So + S) 

< P{<^o{r,),m) exp (P(C:(r?))r) (||uo||?^ 2 iv + || 0 o||^ 2 iv + So + S). 

Thus 


m+2 

(\\dtU\\'^L2ff2(m+l)-2j+3 + ||5iU||^ooj|^2(m+l)-2j+2) + 11 I (.Tt)* 

jf=m+l 

m+2 

+ ^ (^||^P|li2j|^2(m+l)-22+2 + 11 11 ^oo j:^2(m+l)-2i+l ) 


(4.92) 


m+2 

+ ^ (\\d{^\\‘L^H^(m+l)-2j + 3 + ll^^lli !x>Jj2(m+l)-2j + 2 ) + ii8r"«ii(*H. 


< P(^o(^)+(??))exp (P((£(T/))r) (||uo||^ 2 iv + ll^'oll+iv +So + S)- 


Thus, in order to derive the second and third assertions of Pm+ij it suffices to prove 
that 


\ \\dtU\\‘^2}j2{m+l)-2j+3 + ||tI+||^2/^2(m+l)-22+2 + \\d{ 6\\‘^^2 fj2(m+l)-2j+3 

j=0 


+ {\\diu\\l^j^2(m+l)-2j+2 + \\d{p\\ cxD 2 ( m +1) — 2 j+1 + IK«II ^oo^2(m + l) —2j+2 ^ 


n(4.93) 


j=0 


P((£o (??)+(??)) exp (P(£(7?))T) (||uo||+jv + ll^'oll+iv +So + S)- 


In order to prove this estimate, we will use the elliptic regularity of Proposition 13.51 
with k = 2N and iteration argument. As the hrst step, we need the estimates for the 
forcing terms. Combining (I4.76P with the estimates ()4.62p and ()4.63p of Lemma 14.71 
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implies that 

m+l 


E 11'"' 


LJ||2 


l2 ^2m-2j+l 


II 7?3,i ||2 I II p4j ||2 

llL2j|^2m-2j+i + ll-e Ilf 


l2jj2m-2j+3/2 


i=i 


II/?5,i I|2 II 7?l>i ||2 _|_||p3ii||2 

11-^ Il/f2j^2m—2j+3/2 “T ||-f^ \\]^oo J-f2m—2j \ ||-t \\f 

_l_||/?4j||2 _|_||7?5 j||2 

< ||-^ \\J^oo jj2m — 2j + l/2 < ||-^ Hr 


l^co}j2m-2j “t" 11-^ llL°°ff2m-2j 
2 

1^00 Jj2m — 2j+ 1/2 


<p{m){^+ E(ii ^ 11 L2 j|^2m-2£+3 + \\dlO\\ Jj2 J^2m—2f.^3 
^ £=0 

+ ^ (^||9fM|||ooJ|/2m-2<+2 + \\dl9\\\^^2-m-2t+2 + \\dtP\\\'^H'^ 


(4.94) 


!7ti-2.^+2 


£=0 


+ 118,'pII" 


J^oo j^2m — 2^+l 


< 


P(eo(r?),i^(r?))exp {P{^{ii))T) (||no||^2iv + || 0 o||^ 2 iv + 5 o + ^), 


The estimates of (j4.90p , (I4.76P as well as ()4.5ip , ()4.53l) of Lemma 14.61 allow us to 
deduce that 


7fL>)”u|||ooJ^2 + \\dtD'/^u\\\2H^ 


< 


DTu - DT^\\\i^^2 + WdtDTu - D 


m+l 


U\ 




I II n"'+i,/i|2 I II n^+i-'iis 

+ U\\j^aalj2 “T IIM 


L2/f3 


(4.95) 


< 


Pim) (\iDr4l~„2 + \\DTnf^,„>) + iid: 


r‘< 




+ II A" 


m+l ||2 


< P((£o(r/),J^(??))exp (P((£(77))r) (||uo|||-2iv + ||6'o|||-2iv + 5 'o + 5 ')- 

Since ^4.741) is satisfied in the strong sense for j = m, for almost t G [0, T], {D/^u, d/^p, 
d/^0) solves elliptic system (13.11) with replaced by — dtD/^u, = 0, F^ 
replaced by — dt{d/^6) and F^, F^ replaced by F^’™', F^’™', respectively. Then, 
we apply Proposition 13.51 with r = 5, then square the resulting estimate and integrate 
over [0,T], to deduce that 

II A”^^iii2^5 + ii9rpii2^4 + ii5r0iii2^5 

< ||Fl’™ - dtDTu\\l2H3 + ||i"'’™ - dt{dT9)\\l2H3 

I II 774,m||2 I II 77'5,m||2 

11“^ llL2i(f7/2 “T \\]^2jj7/2 

< 11^1,m||2^^^ ||(9fFj"u||^2j|^3 + ||T^’™'|li2j|^3 + \\dt{d/^9)1112^3 

I II 77'4,m||2 I II 77i5,m||2 

■r IT IT2JI/7/2 i- IT IT2 j|^7/2 

< F((2:o(r/), 1^(7?)) exp (F(e(r?))r) (||uo||^ 2 iv + ||0o|ll,2iv + + (?), 


(4.96) 


where in the last inequality, we have used ^4.901) . p4.94p and (I4.95p . Similarly, Propo¬ 
sition 13.51 with r = 4 reveals that 

II A”^^iiioc^4 + ii9rpiiioc^3 + iiar0iiioo^4 

< iiFi’™ - 9iFr^iiioc^2 + iif3-- - dt{dr9)\\i^H2 

I II 774,m||2 I II 77’5,m||2 

“•“11-^ \\l^oofj5/2 -r \\^ \\l^oofj5/2 

< F(go(r/),i^(r?))exp(F(e(r?))r) (||uo||^2iv + ||0o|ll^2iv + + (?)• 


(4.97) 
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By iterating to estimate d^u, d^p and dlO for j = 1,... ,m, as well as the above two 
estimates, we have that 

< P{(ao{v),m) exp {P{^{v))T) {WuoWl,^ + \\ 9 o\\l 2 N +do + d). 

Thus, we have that 


E 

i=i 

+ 


^^11 l2fj2(m+l)-2j+3 + \\dtP\\ 1^2 fj2{m+l)-2j+2 + Hn ]^2jj2(m.+ l)-2j+3 


m 

E 

i=i 


11 7^00 ^ 2 (m+l) — 2 j + 2 + II^pII 7 ^ 007 ^ 2 ( 771 + 1 ) — 2 j + l + Ili¥e|l 7^00 77 ' 2 (m+l) — 27 + 2 


/4.98) 


< 


P{(Bo{r]),Si{r]))exp {P{(S{r]))T) (||uo|||- 2 iv + WOoWh^n +5o TS")- 


Then we apply Proposition 13.51 with r = 2(m + 1) + 3 < 2N + 1, square the result 
estimate and integrate over [0, T] to see that 

\W\\ ^2 f{2(m+l)+3 + \\p\\ 12 f{2(m+l)+2 + 11 ^ 11 ^,2 J|^2(m+l)+3 


< 


\F dtU\\j^2ff2(m+l) + l P \\P dtO\\j^2ff2(m+l) + l 


+ l|F 


4||2 


\l2jj;2(m+l)+3/2 


+ ||F®ii 2 


l2jj;2(m+l)+3/2 


< II pl l|2 

^ IT IT 2 j:^ 2 (m+l) + l 


+ II^tII 


I2ff2(m+1) + 1 


+ l|i^' 


3||2 


I ]^2H2(m+l) + l 


+ \m\ 


(4.99) 


I2jj2(m+1) + 1 


+ l|F 


4||2 


I l2ff2{m+l)+3/2 


+ l|i^' 


'5||2 


I l2ff2{m+l)+3/2 


< 


P{<Bo{v),m)^W {P{m)T) {\\uo\\l2N + ll^olln^iv +5o + 5), 


and then again with r = 2(m + 1) + 2 < 2N to see that 

ll'^ll7;^oo77-2('m,+ l) + 2 “1“ ||p||7(^oo77'2(m+l) + l + 11^11 7^00 71^2(771+1)+ 2 
< ||F^ — 9itt||^ooj|^2(m+l) + \\F^ — dt9\\‘^j^^jj2(7n+l) 


+ \\F- 


4||2 


+ l|i^' 


< 


J^co fj2(m+l) + l/2 
1||2 , lln„.l|2 


5||2 


I L<xi fj2(m+l) + l/2 
P3||2 


F 11 ^ 00 ^ 2 ( 771 + 1 ) + 7)^00 77'2(777 + 1) H“ 11-^ 11 77oo77'2(m+l) + 119+11 77oo77'2(7ti+i) 


(4.100) 


I II p4 ||2 ill pS ||2 

~F r 00 Lr2(’m+l)+l/2 ' 


TooJ|^2(m+l) + l/2 “T 11+^ |l2^ooJ|^2(m+l) + l/2 

< P{(Bo{p), ii{rj)) exp {P{(B{p))T) {\\uo\\]^ 2 n + || 6 'o|||- 2 iv +5'o + (?)- 
Thus (j4.93p is obtained by summing (j4.98p - (j4.100p . This completes the proof. □ 


5. Preliminaries for the nonlinear problem 

In order to use linear theory for the problem (11.811 to solve the nonlinear problem 
», we have to define forcing terms F^, F^, F^, F^ to be used in the linear estimates. 
Given u, 9, rj, we define 

F^{u,9,r}) = dtfi{l + xz)Kd 2 ,u — u-V a/u and F'^{u,9,ri) = p^yV, 

F^{u,9,r]) = dtfi{l + X3)Kd39 - u ■Vj^/9 and F^{u,9,r]) = - {py'l, 


(5.1) 
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where , JV^ K are determined as before by 77 . Then we define the quantities 
^n{u-,0) and ^n{u) as 

N 

^n{u, 0) = {\\9tU\\\2jj2N-2j+l + \\d{u\\‘j^aojj2N-2j 

J=0 ^ (5.2) 

+ \\dtG\\\2H2N-2j+l + \\d{0\\‘j^^^2N-2j'j , 

and 


N 


^ y\\d{u\\‘^j^2fj2N-2j+l + \\dt'^\\‘L°°H^N-2j 

j=0 


(5.3) 


5.1. Initial data estimates. Since r/ is unknown for the full nonlinear problem, and 
its evolution is coupled to that of u, p and 6, we must reconstruct the initial data to 
contain this coupling, only with uq, Oq and rjQ. Here we will define some quantities 
which have minor difference from [6]. 

^0 ■= ||■wo||^2Jv + ||0o||^2iv + ||?7 o||^2jv+i/2) (5-4) 

and 

N N-1 N 

^o{u,P,0) := ^ ||5jM(0)||^2iV-2j + ^ \\d{p{0)\\‘jj2N-2j-l +y~] ||5 j0(O)||^2JV-2j . (5.5) 
j=0 j=0 j=0 

For j = 0,... ,N — 1, 

diiF\u,p, 9),F^u,p, e),F\u,p, 9), F^iu,p, 9)) 

■= ^ ^||9fF^(0)||^2JV-27-2 + \\dfF^{0)\\‘^jj2N-2t-2 + 11 <9^ (0) 11 ^2JV-27-3/2 (5.6) 

£=0 

+ \\dtF^{0)\\‘^^2N-2e-3/2 

and for j = 1 


^o(^) •— ll%ll^2iV+l/2, 




^oiv) ■— ll%ll^2iV+l/2 + ||'9t^(0)||^2JV-27+3/2 ■ 

£=1 

(Sq(u,P, 0) := ||tto||^2iV + ||0o||^2iV, 


(5.7) 

(5.8) 

(5.9) 


and for j = 1 






(BI{u,p,9) ^ ||9fM(0)||^2iv-27 + ^ \\dfp\\‘lj2N-2t-i + ^ \\df9{0)\\‘jj2N-2e- (5.10) 

e=o £=o 1=0 

The following lemma is a minor modification of Lemma 5.2 in [6], so we omit the 
details of proof. 

Lemma 5.1. For j = 0,..., N , 

\\d{u{0) - Dlu{0)\\‘lj2N-2j < Pj{(Sl{r]),(Bl{u,p,9)) (5.11) 

and 

\\di{9{0)V^,y3p) - Di{9{0)V^,y3,o)\\l2N-2, < P,{(Bi{v),(Bi{u,p,9)) (5.12) 
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for Pj{-, •) a polynomial such that Pj{0, 0) = 0. 

For F^{u,9,r]), F^{u,6,r]), F‘^{u,9,r]) and F^{u,9,rj) defined by (15.11) and j = 
0 ,..., — 1, we have that 

di{F\u,p,9),F^iu,p,9),F\u,p,9),F^{u,p,9))<Pj{4-^\r,),<Bl{u,p,9)) (5.13) 

for Pj{-, •) a polynomial such that Pj{0, 0) = 0. 

For j = 1,..., N — 1, let F^’PO), F^’PO), F'^’PO) and F^’PO) are determined by 
dMlD, (irail and dSI]). Then 

\\^^’H^)\\'H2N-2j-2 + ||F^’-^'(0)||^2iV-2j-2 

+ \\F^’H0)\\W2,-^/2 + \\F^’H0)\\l2.-2,-^/2 (5.14) 

<Pfi<Bl+\p),4iu,p,9)) 

for Pj{-,-) a polynomial such that Pj{0, 0) = 0. 

For j = 1,... ,iV- 1, 

2 

< Pj{^l{r]),<Bl{u,p,9)) (5.15) 

jj2N-2j+i/2 

for Pj{-,-) a polynomial such that Pj{0, 0) = 0. Also, 

11^0 ■ '^||^2iV-l/2(j]) < ||■Wo||^2JV (l+||%||^2iV+l/2). (5.16) 

This lemma allows us to construct all of the initial data 5^tt(0), (9^0(0), 9^f?(0) for 
j = 0,..., and dip{0) ioi j = 0,N — 1. 

Assume that Jq < oo. As before, we will iteratively construct the initial data, but 
this time we will use Lemma [5Tl We define dtrj{0) = uo-Pfi), where uq G Lf^^“^/^(S), 
and ^ is determined by r]Q. (|5.16p implies that ||9fi?(0)||^2iv-i/2 ^ Pi^o) for a 
polynomial F(-) such that F(0) = 0, and hence that <^q{u,p,9) + ~ P{<^o)- 

Then (15.131) with j = 0 implies that 

do{F^{u,P,9),F^{u,p,9),F‘^{u,p,9),F^{u,p,9)) < Fo(<£o(??), <^o{u,P,0)) < P{^o) 

(5.17) 

for a polynomial F(-) such that F(0) = 0. Note that in these estimates and in the 
estimates below, the polynomial F(-) of S'q are allowed to change from line to line, 
but they always satisfy F(0) = 0. 

In this paragraph, we will give the iterative definition of d{p{0), +'u(0), cfi+^9{0) 

and for 0 < j < N — 2. Now suppose that dfu{0), df9{0) are known for 

£ = 0,... ,j, dfr]{0) is known for ^ = 0 ,..., j + 1, dfp(0) is known for £ = 0 ,..., j — 1 
(with the exception for p(0) when j = 0) and 

el(u,p,9) + £^■^^(7?) 

+ Si(F\u,p,9),F^iu,p,9),F\u,p,9),F^(u,p,9)) (5.18) 

<Fi^o)- 

And according to (I5.14h and (15.lip , we know that 

\\F>t'u{0)\\‘jj2N-2j + ||F^’'^(0)||^2JV-2i-2 + ||F^’-^'(0)||^2JV-2j-2 
+ l|i""’''(0)||^2^_2,_3/2 + \\F^’H0)fH2N-2,-3/2 

< Pi^o). 


^p)af^(o)-a^-'u(o) 
£=0 ^ 


(5.19) 
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By virtue of estimates (I4.68P 

+ ||f(Fi’^(0),I)^\r(0))||^,,_3._5/. (5.20) 

< Pi^o) 

This allows us to define ^p(O) as the solution to (I3.18P with /^, /^, replaced by 
f^, f^, f. The Proposition 2.15 in [13] with k = 2N and r = 2N — 2j — l implies that 

\\&ipm\mN-2,-.<P{<^0). (5.21) 

Now we define 

9^+i0(O) = £°2^9^-0 (o),f3J(O)) G H^N-2j-2_ (-5 22) 

Then according to (I5.18P and (15.191) . we have that 

||a^'+'0(O)||^,^_3,_. <P(<fo). (5.23) 

Now the estimates p4.66p . ^5.181) and (I5.19P allow us to dehned 

Zl^+^(0) := (f1’^(0) + di{eiO)V^,y3,o),DiuiO),dlp{0)) G (5.24) 

and then according to (15.111) . we have 

\\4^\m\]i2N-2i-2 <P{<^0)- (5.25) 

Now the estimates p5.16p . ^5.181) and p5.25p allow us to dehne 

= E (■' t ■ a(’"‘"'»(0), 

£=0 ^ '' 

and imply the estimate 

\\Pi^^mfH2M-.,-,/.<P{^o). (5.26) 

Thus, (|5.18l) together with (|5.21l) - (|5.26l) imply that 

and then p5.13p implies that 

do^^{F^iu,p,9),F^{u,p,9),F'^{u,p,e),F^{u,p,9)) < F{S’o). 

Hence that we can deduce the estimate 
<Bi+\u,p,9) + 4+\y) 

+ ^i+\F\u,p,9),F^{u,p,9),F\u,p,9),FHu,p,9)) 

< P{<Po)- 
For j = N — 2, we have 

<^o~^{u,P,9) + <^o{y) 

+ S^-\F\u,p,9),F\u,p,9),F\u,p,9),F^iu,p,9)) (5.27) 

< P{<^o)- 

Then, we only need to dehne d^~^p{0), d[^9{0) and d^u{0). Like the construction 
after Lemma 14.81 we need the compatibility conditions on uq and r]Q. Now we have 
constructed d{p{0) for j = 0,...,N — 2, d{u{0), d{9{0), F^’^{0), F^’fO), F^’fO), 
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^^’■^ (0) for j = 0,..., — 1, and dirj{0) for j = 0,..., N. We say that uq and rjo 

satisfy the Wth order compatibility conditions if 

• (-^t“(0)) = 0 in hi, 

< Dlu{0)=0 onSfo, (5.28) 

Ho =0 onS, 

for j = 0, — 1, where Ho is the projection defined as in (I4.14p and Dt be 

the operator defined by (14.101) . Note that if no and % satisfy (j5.28p . then the j-th 
compatibility condition p4.71l) is satisfied for j = 0,...,A^—1. Then the construction 
of d^~^p{0) is the same as [6] using the compatibility condition (|5.28l) and the elliptic 
theory of .e/- Poisson equations (13.181) derived by Y. Guo and I. Tice in P and L. 
Wu in m- And 

l|5f-Vo)||2,i <P(<fo). (5.29) 

Then we set d[^6{0) = 6{0), (0)) e due to (14.671) and (I5.14p . and 

set D^uiO) = + af-^( 0 V. 4 ,y 3 ,o), A^"^^^(0),5f-^p(0)) € due to 

(|4.66p and Lemma 15.11 And D^u{0) G ‘3^ (0) is guaranteed by the construction of 
d^~^p{0). As before, we have 

||9f u(0)||^o + ||5f 0(0)11^0 < P{^o). (5.30) 

This completes the construction of initial data. Then summing the estimates (15.271) . 
(I5.29P and (I5.30p . we directly have the following proposition. 

Proposition 5.2. Suppose that uq, 9q and po satisfy S’q < oo. Let the initial data 
diu{0), di 0{0), dip{0) for j = 0,... ,N and dip{0) for j = 0,... ,N — 1 be given as 
above. Then 

S’o < <Bo{u,p,e) + iBo{p) < P{S’o). 

Here ^o{p) = {p), which is defined in (14.501) . 

5.2. Transport equation. Here we consider the equation 

dtp + uidip + U 2 d 2 P = U 3 on S, 

7/(0) = po- 

The local well-posedness of (j5.32l) has been proved by L. Wu, which is the Theorem 
2.17 in |13] . The idea of his proof is similar to the proof of Theorem 5.4 in [6]. In 
m, L. Wu has proved in Lemma 2.18, that the difference of p and pQ in a small time 
period is also small. 


(5.31) 


(5.32) 


5.3. Forcing estimates. In the next section for the estimates of full nonlinear prob¬ 
lem, we need some forcing quantities. Besides and jJo which have been defined in 
(|4.49p . we define the following quantities 


Af-l 


:\N ttS ||2 


^ ^ {\\9iF' \\x^2H2N-2j-l + \\dtP \\]^2H2N-2j-l) + || 5 ^ P \\H° + \Pt ^ \\L'^H° 

3=0 
N 

+ E 

3=0 


?5ll2 


llL°°iT2iV-2j-l/2('5^) + 11 ^ A" \\j^^jj-2N-2j-l/2(^Y:) I I 
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Af-1 


T-L : — ^ ^||c^'F^||^ 2 j|^ 2 JV- 2 j -1 + \\diF'^\\\2H2N-2j-l 
i=o 


Af-1 

+ E(||^' 

j=0 


i2^2JV-2j-l/2(s) 


+ WdiF^ 


li2^2iV-2j-l/2(s) i ; 


The following theorem is similar to Theorem 2.21 in m with obvious modification. 


Theorem 5.3. The forcing terms satisfy the estimates 

d < P{m) + P{^N{u,e)), (5.33) 

< P{<^o), (5.34) 

p < P{m) + P{^N{u,e)), (5.35) 

H < T {P{^i{r^)) + P{SiN{u,e))) ■ (5.36) 

Proof. The proof of this theorem is the same as the proof of Theorem 2.21 in m , so 
we omit the details here. □ 


6 . Local well-posedness for the nonlinear problem 

6.1. Construction of approximate solutions. In order to solve the (11.41) . we will 
construct a sequence of approximate solutions (u"*,p™', 0 ”*, 77 "*), then take the limit 
m —>• 00 . First, we construct an initial pair {u^,6^,r]^) as a start point, then we 
iteratively dehne all sequences (tt™',^”*, 0 ™, r/™) for m > 1 . 

Suppose that the initial data (no, 0 o 5 0 o) has given. According to the Lemma A.5 
in there exist and 0 ° defined in 11 x [ 0 ,oo) with d{u^{0) = d{u{0), d{9^{0) = 
5^0(0), for j = 0,..., N, satisfying 

SiN{u^,e^)<P{<^o)- ( 6 . 1 ) 

Then we consider the equation ()3.18p with u replaced by u^. From the Theorem 
2.17 in [13], the hypothesis of which is satisfied by (I5.31|) and (j6.ip . there exists a rj^ 
defined in H x [0, Tq), which satisfies 5^ 77 ^( 0 ) = d^r]{0) for j = 0,..., N as well as 

< Pi^o). 

Then for any integer ttt- > 1, we formally define the sequence (u'^, , 9^, rj'^) on 

the time interval [ 0 , T^) as the solutions of system 

' dtu^ - 

= + X3)K^-^d3U^-^ - inll, 

div_j/m-i u™'= 0 inll, 

< dt9^ - A^^- 10 "* = + X3)K^-^d39^-^ - • V^^- 10 ™-^ inll,( 6 . 2 ) 


S^m-i{p^,u^)PK^-^ = onS, 

+ = -\jT^-^\ onS, 

= 0, 0”" = 0 on Sb, 

and 

= onS, (6.3) 


where A™ ^ are determined in terms of 77 "* ^ and is in terms 

of 77 ”*, with the initial data ( 7 X^( 0 ), 0 ^( 0 ), 77 ”^( 0 )) = (tto, 0 o,ho). 
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In the following, we will prove that these sequences can be defined for any integer 
m > 1 and the existence time Tm does not shrink to 0 as m —)• oo. The following 
theorem is a modihed version of Theorem 2.24 in m, which improves the estimate 
()4.75l) using the energy structure and elliptic estimates. 

Theorem 6.1. Suppose J(0) > <5 > 0. Assume that the initial data (uo,6*o,??o) satisfy 
(§0 < oo and dfU(O), 8^0(0), d^ri{0), for j = 0,...,N, are given as above from the 
Proposition \5.S[ Then there exists a positive constant ^ < oo and 0 < T < 1 depend¬ 
ing on S’q, such that if 0 < T < T, then there exists a sequence 

(when m = 0, the sequence should be considered as {u^, satisfying the iteration 

equation (|6.2I) within the time interval [0, T) and the following properties: 

1. The iteration sequence satisfies 

An{u'^,0'^)+A{t]'^) < (6.4) 

for any integer m>0, where the temporal norm is taken with respect to [0,T). 


2. J'^{t) > 5/2 with 0 < t < T, for any integer m > 0. 

Proof. In this proof, we will follow the path of proof of Theorem 2.24 in |13] . We 
will use an inhnite induction to prove this theorem. Let us denote the above two 
assertions as statement Pm- 

Step 1. Pq case. The only modification here is that the construction of and 0^ 
reveals that A.]y{u^,0^) < P{^o). Then the rest proof of this case is the same as the 
proof of Theorem 2.24 in m- Hence, Pq holds. That is A.n{u^, 0^)-\-R{r]^) < J" with 
the temporal norm taken with respect to [0, T) and J^{t) > 5/2 for 0 < t < T. 

In the following, we suppose that Pm-i holds for m > 1. Then we will prove that 
Pm also holds. 

Step 2. Pm case: energy estimates of 0™ and u'^. By Theorem 14.91 the pair 
p"^, 0'^) satisfies the equation 




- inH, 

{D/^u'^) = 0 inH, 

< dt{d/’0^) - 0^) = inH, 

S^m-i p”*, on S, 

V^m-i{d/^0^) ■ = F^’^ onS, 

^D/^u^ = 0, 5f0™ = O on Eft, 


(6.5) 


in the weak sense, where F^’^, F^’^, F^’^ and F^’^ are given in terms of p™, 
0"*, and 0™'“^, Then for any test function </> £ where 

(J^)™”^ is the space with r] replaced by rf^~^, the following holds 
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Therefore, when taking the test function (j) = 9^, we have the energy structure 

^ JQ. Jo JQ, Jo JTl 

= 1 f J™"^(O)|5f0”^(O)|2 + i f f (6.6) 

2 Jo 2 Jo Jo 

+ r / + f ! F^’^d^e^. 

Jo Jo Jo Je 

By induction hypothesis, (|5.31l) . trace theory and Cauchy inequality, we have 

Wdfe^Wl^HO + ||9f 

< sup (i / J— r / r / 

0 <t<T\^Jo Jo Jo Jo Je / 

2 Jo 2 Jo Jo Jo Jo 

+ [ [ F^’^dtO'^ (6-7) 

Jo Je 


< P(<fo) + T^||5f + Vr||F3-^||i,^o 


for a polynomial P{0) = 0. Taking T < min{l/4,1/(165°^)} and absorbing the extra 
terms on the right-hand side into left-hand side imply 


l|5f 0”^llioc^o + Wdfo^wl,^, < P{So) + Vt\\f^’^\\1,^, + VT\\F^^^\\l^^_,„^^y 

( 6 . 8 ) 

By induction hypothesis, we have 


I II2 

T Wl^ho 


' N-1 


< 


< 


P(i^(r7™-i)) mu^WhH. + \\di9^\\hH. I +F 
\j=o 

P(<fo + 3f)+F, 


p5,Af||2 

^ IIl°°_h'-i/2(s) 

\\d;u”'\\l«„, + \\din\i^HA 

<P{£'o + ^)+T. 

And, the energy estimates about vF is the same as the proof of of Theorem 2.24 in 
m- Therefore, we have 

Wdfu^Wl^Hi + l|5f 0™||i2^i < P{^o) + ^/TP(J’o + 2f) + VfF. (6.9) 


P{Si{r] 


m—1 
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Step 3. Pm case: elliptic estimates for 0”*, u^. For 0 < n < N — 1, the n-th order 
heat equation is 

(dtid^e^) - inn, 




= F 


5,n 


= 0 


on S, 
on Sft. 


( 6 . 10 ) 


The elliptic estimate in the proof of Lemma 13.31 reveals that 

\\rYnnm\\2 ^ ||773,n||2 , ||on+lflm||2 ,|| 775,71112 

~ ll-C \\ 12 }J2N-271-1 ^ W^t " 11 L2/^2iV-2n-l “T 11 ^ ||2^2j|/2iV-2ri-l/2 • 

( 6 . 11 ) 

As what we did before, 


I II2 

\-^ W 1^2 }J2N-2n-l 


' N-2 


< 


TP{K{ri^-^)) Y. \\dien\UH2.-2,-i + ||9^u™||i^^2.-27- 
\i=o 

<rp(7fo + ^)+H. 

I Z?5?^l|2 

\-^ \\]^2u2N-2n-l 


+ n 


'N-2 


< 


rp(7!(, 

\i=o 

<TP{S’o + 3f)+n. 


+ n 


But for the term 6^1112H2N-2r7-i, we estimate backward from — 1 to 0. First, 
when n = N — 1, this is the case of energy estimate of 9^. Then we iteratively use 
the elliptic estimates (16.lip from n = N — 2 to n = 0 to obtain all the control of 

II an+lam||2 

l|C'i O \\j^2fj2N-2n-l- 

And the elliptic estimate for n™ is the same as the proof of of Theorem 2.24 in 
m- Thereore, we have that 

Af-l 


E (ii«i 


n^.m\\2 


U 


I ^2 J^2iV—2n+l 


+ ii^r^ 


in^m||2 

l2}j2N-2n+l ) 


„.o (6.12) 

< P(<?o) + VfP{l +g, + ^) + VfP + H. 

Step 4. Pm case; synthesis of estimates for vF and 0™. Combining (|6.8I) . (16.121) and 
Lemma 2.19 in m, we deduce that 

SiNiu"^, 9^) < P(^o) + VtP {<^0 + ^) + VtF + n. (6.13) 

Then by the induction hypothesis and the forcing estimates of Lemma 15.31 we have 
that 

F < Pi^iv^-^)) + P{^N{n^-\9^-^)) < P{^), 
n<T {P{R{rj^-^)) + P{iiN{u^-\9^-^))) < TP{3f). 

Hence we obtain the estimate 


9^) < C (p{^o) + VtP{,^o + J^)) 


(6.14) 


for some universal constant C > 0. Taking > 2CP{S‘q) and then taking T sufficient 
small which depends on we can achieve that 0™) < 2C'P((fo) < ^■ 







46 


YUNRUIZHENG 


Step 5. Prn case: estimate for r/™ and These estimates are exactly the same 

as the proof of of Theorem 2.24 in m- So we omit the details here. 

Thus, we can take = P{Sq) for some polynomial P{-) and T small enough 
depending on to deduce that 


^n{u^,0^)<3^ (6.15) 

and 

J^{t)>d/2 forte[0,r]. (6.16) 

Hence Pm holds. By induction, P„ holds for any integer n > 0. □ 

Theorem 6.2. Assume the same conditions as Theorem \6.1[ Then 

R{u^,p^,9^)+K{vn < P{^o) (6.17) 

for a polynomial P{-) satisfying P{0) = 0. 

Proof From the estimates (j4.75p . (I5.3ip . Lemma 15.31 as well as Theorem 2.17 in |13] . 
we directly have that 

R{u^,p^,e^) + Air^^) < p(^o) + p{iiN{u^,en + Hhn)- 

Then, applying the Theorem 16.11 we have that 

A{u^,p^,e^)+A{ri^)<P{<^o)- 


□ 


6.2. Contraction. According to Theorem 16.21 we may extract weakly converging 
subsequences from {(u"*,p”^, 0™, Unfortunately, the original sequence {(u”^,p™', 0'", ? 7 ”^)}“^q 

could not be guaranteed to converge to the same limit. In order to obtain the de¬ 
sired solution to (11.41) by passing to the limit in (16.21) and (|6.3I) , we need to study its 
contraction in some norm. 

For T > 0, we define the norms 

^{v,q,Q-,T) = ||u|||ooj:^2 + ||u|||2j|^3 + WdtvWlooffO + \\dtV\\‘l2ffl + llgll^ooj^l + \\q\\‘l2jj2 

+ I|0|Il°°//2 + II®IIl2j:^3 + \\dtQ\\\oafjO + ||'9t0||^2^1 (6.18) 

m{C,T) = ||C|liooH5/2 + + 115^(1112^1/2, 

where the norm PPH^ is LP([0, T]; i7^(H)) in ill, and is LP([0, T]; i7^(S)) in 911. 

The next theorem is not only used to prove the contraction of approximate solu¬ 
tions, but also used to verify the uniqueness of solutions to (|1.4p . To avoid confusion 
with {(u™,p™', 0™, 77 ”^)}, we refer to velocities as , pressures as q^, temperatures 

as 0-' , , and surface functions as for j = 1,2. 

Theorem 6.3. For j = 1,2, suppose that , q^, 0-^, and satisfy the 

initial data d^v^iO) = clfu^(O), (9j'0^(O) = for k = 0,1, g^(0) = g^(0) and 
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C^(0) = C^(0), and that the following system holds: 


' dtv^ - = dtC^{l + x:,)K^d^w^ 


- ■ V 

inPt, 

O 

II 

inPL, 

dtQ^ - A^,0^' = dtC^{l + X3)K^d3d^ - • V^jd^ 

inPt, 

fp 

\ 

!o. 

Co. 

Co. 

X 

Co. 

II 

Co. 

Co. 

onTi, 

v^i0^' • + 0^' \pt^\ = -\jp^\ 

onTi, 

O 

II 

© 

o' 

II 

onTih, 

II 

IT 

on's, 


where are determined by . Assume that A.{v^,q^ ,Q^), A.{w^ and 

A{C^) are bounded by 2^. 

Then there exists 0 < Ti < 1 such that for any 0 < T < Ti, then we have 

q^ - 0^ - 0^; T) < - w^,0, - d^;T), (6.20) 

9Jt(C^ - C^; T) < - w^,0, - d^; T). (6.21) 

Proof. This proof follows the path of Theorem 6.2 in [6]. First, we define v = 
w = — w'^, 0 = 0^ — 0^, d = d^ — d'^, q = q^ — q^. 

Step 1. Energy evolution for differences. Like the proof of Theorem 6.2 in [6], we 
can derive the PDE satisfied by v, q and 0: 


' dtv + div^i S^i{q,v) - 0V^iy3 = div_j/i(B(^i_^ 2 )ti^) + 

inD, 

div^i V = H'^ 

inD, 

dtQ - A^10 = div^i(V(^i_^2)02) + 

inD, 

S^i {q, 

onS, (6.22) 

V^i0 • + 0 |=yKi| =-V(^i_^2)02 . 

on S, 

II 

O 

© 

II 

o 

onSfe, 

^ v{t = 0) = 0, &{t = 0) = 0, 



and the PDE satisfied by dtv, dtq, dt& from taking temporal derivative for the above 


system: 

' dt{dtv) + div^i S^i{dtq,dtv) - dt{QV^lyl) 

= div^i(Dgj(^i_^2)r;^) + inD, 

diYdtv = infl, 

dtidtO) - A^idtQ = div^i(V(gj^i_a^^2)0^) + inD, 

S^i {dtq, dtv)J^^ = on S, 

V^iai0-^^ + at0|^i| =-Va,(^i_^ 2 ) 02 .^i + #5 onS, 
dtv = 0, dt& = 0 onSfe, 

^ dtv{t = 0) = 0, dtQ{t = 0) = 0, 


(6.23) 
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where and have been given by Y. Guo and I. Tice in [B], 

= e^V^i_^2yl + eW^2{yl - yl) + div^i_^2(B^2U^) - V^i_^2q‘^ 

+ + X3)K\dsW^ - d3W^) + {dtC^ - dtC^){l + X3)K^d3W^ 

+ dtC^(l + X3){K^ — K^)d3W^ — {w^ — tv^) ■ — v? ■ Vj^i{w^ — w^) 

- ■ V^1_^2W^, 

= div^i_^2(V^2 02) + + X3)K^{d3^^ - 

+ - dtC^){l + X3)K^d3^^ + dtC\K^ - K^)d3W^ - {w^ - w^) • 

- • V^i(79^ - _ yj^ . V^l_^2^\ 

= -V^202 . (^1 - - 02 (|^i| - |^2|^ , 

= dtH^ + divg^^i(B^i_^2u2) + div^i (B^i_^25tu2) + div^^^i (B^iu) 

+ div^i(Ba^^iu) - Vq^^iq, 

= dtH^ + divg^^i(V(^i_^2)0^) + div^i(V(^i_^2)(9t0^) + divg^^i V^i0 
+ div^i Vg^^l0, 

= dtH^ - V(^i_^2)9402 . - V(^i_^2)02 . - V^i0 • 

Then we can deduce the equations 

i f \dtv\'^J^{t) + ^ j f \B^idtv\^J^ 

J t/ 0 'J 

= l f [ \dtv\\dtJ^K^)J^+ f f dt{evW3)-dtvJ^ 

^ Jo Jn Jo Jn 

+ [ [ j\H^-dtV + H^dtq) 

Jo Jn 

-If [ - [ [ H^-dtv, (6.24) 

^ Jo Jn Jo Jt. 

\ [ \dte\^j\t)+ f [ |V^i5i0|Vi+ f [ |5t0|2|^i| 

= l f f \dte\\dtJ^K^)J^ + f f • dtQ 

^ Jo Jn Jo Jn 

- t [ ■/'Va,^i_a,^202.V^i5i0+ r [ H^-dtG. 

Jo Jn Jo Jt, 

Step 2. Estimates for the forcing terms. Now we need to estimate the forcing terms 
that appear on the right-hand sides of ()6.24p . Throughout this section, P{-) is written 
as a polynomial such that T’(O) = 0, which allows to be changed from line to line. 
The estimates for ||^^||o, ||.ff2||Q^ ||5t.ff2||p^ ||^^||_i/ 2 ) \\H%, \\H%+i, \m\ r+1/2, 
II diVj^i(B(^i_^2)u2)||^ and ||B(j^i_^2)u2^^||^_,_]^y2 have been done by Guo and Tice 
in [6]. So we can directly using them only after replacing e by By the same 
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method, we can also deduce that 

ll-^^llo ^ -P(V^)(||0||2 + IIC^ - C^ll3/2 + - dtQ^\\l/2 + \\dtC^ - '9t^C^lll/2,g 25\ 

+ ||«i‘ - oj'llo + lla.to' - a.totlo + l|tf‘ - tftli + 

\\H^\\-in ^ Av^)(IIC‘ - C"lli/2 + lia<C' - 8 ,C|Ii/ 2 + lieib), 

and for r = 0,1, 

< P{^){\\C^ - C2||,+i/2 + - dtC\-l/2 

+ llw^ - w'^Wr + ||?9^ - I?^||r+l), 

ll^^'^ll.+l/2<^’(^/^)IIC'-C'L+3/2, 

II div^i(V^i_^.02)||, < P(V^)||C' - C'||2.+3/2, 

||V^i_^ 202 . ^l||,+i/2 < n^^)IIC' - C'llr.+3/2. 

Step 3. Energy estimates of dtv and 5*0. First, owing to the assumption and 
Sobolev embeddings, we obtain that 

||ji||Loo + ||Ki||ioc <1 + P(v^) and IlStJ^IU- <P(v^). (6.31) 


(6.26) 

(6.27) 

(6.28) 

(6.29) 

(6.30) 


The bounds of (|6.31l) reveals that 

Iff \dtQ\\dtJ^K^)J^<P{V^)l f [ \dtQ\\j\ (6.32) 

^ Jo Jn ^ Jo Jo, 

In addition, estimates (I6.25D . (16.261) together with trace theory and the Poincare 
inequality reveals that 


• 5*0 - 

< f [ II^'IIl-(|I^'IIt- 

+ /V'|I-1/2||5*0||i/2 
Jo 

< [ P{^)^/Z, 

Jo 




. V^i5*0 - f [ H^- 5*0 

JO JE 

||.ff^||o||5*0||o + ||Vaj^i_aj^20^||o||V^i5*0||o^ 


(6.33) 


where we have written 


2: = lie - eiii/2 + ii5*e - 5*eii?/2 + - d^,cYp2 

+ ||u;^ - w'^Wl + \\dtw^ - 5iu;2||2 + ||*9i - *?2||2 ^ _ g^^2||2 ^g_34) 

+ ieili + lkll? +110111- 


Combining (|6.32l) . (|6.33l) . (|6.24l) . Poincare inequality of Lemma A. 14 in [6] and Lemma 
2.9 in [13] and utilizing Cauchy inequality to absorb ||5*0||i into left, yield that 

^|8,e|Vi(i) + ij(‘||8,0|i; 

<F(\/5)i / |8,e|V‘(()+ j P{sf^)Z 
^ Jn Jo 


(6.35) 
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Then Gronwall’s lemma and Lemma 2.9 in m imply that 

WdtOWl^jjo + 115*011^2^1 < exp{P(^/^)T} r P{^)Z. (6.36) 

Jo 

Then energy estimates for dtv are likely the same as what Guo and Tice did in [6], 
so we omit the details. The energy estimates for dtV and 5*0 allow us to deduce that 

||5*u||^oo*^0 + ||5 *u||^ 2*^1 + ||5*0||^oo*^O + ||5*0||^2*yl 


< exp{P{V^)T} 


J='(V5)||9|li.„o+C||8,C‘-8iC 

+ - i>l(X oc^l/2 + ||u||*;,oo^ 

^ J=0 

L2m/2 + 


i2*^-i/2 + [ P{yf^)Z 
Jo 


(6.37) 


■j=o 


where the temporal norm of L°° and are computed over [0, T]. 

Step 4. Elliptic estimates for v, q and 0. For r = 0,1, we combine Proposition 
()3.5p with estimates (j6.27p - (j6.30p as well as the bounds of ||77^||r, ||77^||r+i, ||77^||r+i/2 
II div_(/i(B(^i__j/ 2 )U^)||r, ||B(_j/i_^ 2 )U^c/L^||r_|_i /2 done in the proof of Theorem 6.2 in 
[6] to deduce that 

lbllr+2 + lkllr+1 + ll®llr+2 

^ Cim)(^\\dtv\\r + ||div^i(D(^i_^2)*;^)||2 + ||iL^||2 + ||i7^||^+i + ||5*0||2 

+ ll-f^^llr + II div^i(V^i_^20^)||^ + ||D(^i_^2)*;^^^||^_,_j^^ 2 + ll-f^'^llr+ 1/2 

+ l|V^i_^20^ • ^^llr+ 1/2 + ll-f^^llr+ 1 / 2 ^ 3g^ 

~ C{r]o) ^||5*u||^ + ||5*0||^ + IIC^ — C^llr+ 1/2 
+ P(v^)(||C^ - C^II^+ 3/2 + ||5*C^ - 5*C^||2_^/2 
+ 11^^^ - W^^llr+l + - ^9^llr+l)) • 


Then we take supremum in time over [0, T], when r = 0, to deduce 

ll^lli°°H'2 + Ikllioo/*! + ||0|||oo*^2 

~ Civo) ^l|5*u|||ooj*0 + ||5*0 ||^c3o*^O + IIC^ — C^II^OO*1*1/2 
+ P(\/^)(||C^ — C^llioo**3/2 + ||5*C^ — 5 *C^||^oo**-1/2 
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Then we integrate over [0, T] when r = 1 to find 

+ lkll + II®IIl2H3 

~ Civo) (^\dtV\\‘j^2Hl + 1119*011^2^^1 + IIC^ — C^||^2j:*3/2 
+ p(v^)(iic' - c'iii2^5/2 + - 5*c'iii2^i/2 

+ ||W^ — w‘^\\j^2H2 + ||l9^ — ^?^lli2/*2)^ • 

Step 5. Estimates of and contraction. After making preparations in the 

above steps, we can derive the contraction results. Since this step follows exactly the 
same manner as the proof of Theorem 6.2 in [6], we omit the details here. Hence, we 
get the ()6.20l) and ()6.21l) . □ 

6.3. Proof of Theorem II.IL Now we can combine Theorem 16.21 and Theorem 16.31 
to produce a unique strong solution to ()1.4p . It is notable that Theorem 11.11 can be 
directly derived from the following theorem, which will be proved in the same manner 
as the proof of Theorem 6.3 in [6]. 

Theorem 6.4. Assume that uq, Oq, r]o satisfy (fo < oo and that the initial data 
dfu{0), etc. are constructed in Section rOl and satisfy the N-th compatibility condi¬ 
tions (|5.28l) . Then there exists 0 < Tq < 1 such that if 0 < T < Tq, then there exists 
a solution {u,p,9,i]) to the problem (II.4p on the time interval [0,T] that achieves the 
initial data and satisfies 

Si{u,p,9)+A{7])<CP{c^o), (6.41) 

for a universal constant C > 0. The solution is unique through functions that achieve 
the initial data. Moreover, rj is such that the mapping defined by (11.21) . is a 

(j 2 N-i cUff^Qmorphism for each t S [0,T]. 

Proof. Step 1. The sequences of approximate solutions. From the assumptions, we 
know that the hypothesis of Theorems 16.11 and 16.21 is satisfied. These two theorems 
allow us to produce a sequence of 6*"*, which achieve the initial 

data, satisfy the systems (|6.2p . and obey the uniform bounds 

sup iA{u^,p^, 9^) + A{rr)) < CP{,^o). (6.42) 

m>l 

The uniform bounds allow us to take weak and weak-* limits, up to the extraction of 
a subsequence: 

d{u^ ^d{u weakly in L‘^{[0,T]-ioi j = 0,..., N, 

dN+l^m ^ ^7V+l^ 

d{u'^ A weakly — * inL°°([0, T]; (Q,)) for j = 0,..., N, 

d{p^ d{p weakly inL^([0, T]; (H)) for j = 0,..., N, 

dipm A Qjp weakly - * in L°°([0, r]; (H)) for j = 0,..., N, 

di9^ -^d{9 weakly T ]-for j = 0,... , N, 

QN+igm df^^9 weakly in 

&i9^ ■^d{9 weakly - * inL°°([0, T]; ^^^-^^(H)) for j = 0,..., A, 
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and 

^ Qj^ weakly inL2([0, T]; for j = 2,..., + 1, 

rj'^ ^rj weakly-*mL°°{[0,T];H^^+^/‘^{J:)), 

dlrj^ A qIt] weakly - * in L°°([0, T]] for j = 1,..., N. 

The collection {v,q,Q,C) achieving the initial data, that is, dlv{0) = &lu{0), 9^0(0) = 
d{6{0), C(0) = <9/^(0) for j = 0,..., and d{q{0) = d{p{0) for j = 0,..., A^ — 1, is 

closed in the above weak topology by Lemma A.4 in [6]. Hence the limit {u,p,6,rj) 
achieves the initial data, since each 0™, r/”^) is in the above collection. 

Step 2. Contraction. For m > 1, we set 

^ ^ pm+2^ ^2 ^ ^m+l^ 0l^0m+2^ 02 ^ Qm+l ^ ^1 ^ ^m+l^ ^2 ^ 

_ j.ym+1^ (^2 _ ^ Then from the construction of initial data, the initial data of 

, q^, , ‘d^, math the hypothesis of Theorem 16.31 Because of (16.2p . (|6.19p 

holds. In addition, (I6.42p holds. Thus, all hypothesis of Theorem 16.31 are satisfied. 
Then 

0 I(n ”*+2 _u^+\pm+2 _pm+l^ Qm +2 _ Qm+1. 

1 1 1 1 (6.43) 

<-0T(n™+^ - - p"*, 0"*+^ - 0™; T), 

^ gm+l _ gm. j.y ^g_ 44 ) 

The bound (I6.43P implies that the sequence {(u”^,p™', 9^)}'^^q is Cauchy in the norm 
Thus 


(6.45) 



inL°° 

([0,r];772(f^)) 

n 


dtu^ dtu 

inL°° 

{[D,T]-H\a)) 

n 

([0,r];77i(fl)), 

p"^ ^ p 

inL°° 

{[D,T]-H\a)) 

n 



inL°° 

{[D,T]-H\a)) 

nL^ 


^ dt9^ ^ dt9 

inL°° 

{[D,T]-H\a)) 

n 

([0,r];77i(fl)), 

. Because of (|6.44p. wf 

deduce that the sequence 


the norm y^9Jl(-;T). Thus, 


rj^ ^ Tj 


dtV^- 

dh" 


dtp 

► d'fp 


mL°° 

mL°° 

inL^ 


([0,r];F5/2(E)), 

([0,r];F3/2(E)), 

([0,r];iLV2(s)), 


(6.46) 


as m —)• oo. 

Step 3. Interpolation and passing to the limit. This section is exactly the same as 
the proof of Theorem 6.3 in[^, which gives the existence of solutions and the estimate 

(IQTD . 

Step 4. Uniqueness and diffemorphism. This section is similar to the proof of 
Theorem 6.3 in^. □ 
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